arxiv: 0802. 0777 
AEI-2008-003 
ITEP-TH-06/08 



Quantum Deformations of the 
One-Dimensional Hubbard Model 

Niklas Beisert 11 and Peter Koroteev 1 " 1 * 

^ Max- Planck- Institut fur Gravitationsphysik 

Albert- Einstein- Institut 
Am Muhlenberg 1, 14476 Potsdam, Germany 

" Institute for Theoretical and Experimental Physics 
B. Cheremushkinskaya 25, Moscow 117259, Russia 

* Institute for Nuclear Research 
Prospekt 60-letiya Oktyabrya la, Moscow 117312, Russia 

nbeisert , koroteev@aei . mpg . de 



Abstract 

The centrally extended superalgebra psu(2|2) x M 3 was shown to play 
an important role for the integrable structures of the one-dimensional 
Hubbard model and of the planar AdS / CFT correspondence. Here we 
consider its quantum deformation U g (psu(2|2) x IR 3 ) and derive the fun- 
damental R-matrix. From the latter we deduce an integrable spin chain 
Hamiltonian with three independent parameters and the correspond- 
ing Bethe equations to describe the spectrum on periodic chains. We 
relate our Hamiltonian to a two-parametric Hamiltonian proposed by 
Alcaraz and Bariev which can be considered a quantum deformation of 
the one-dimensional Hubbard model. 



1 Introduction and Overview 



Finding the spectrum of a quantum mechanical model is an intricate problem. Indeed, for 
generic models there is no complete analytic solution to the spectrum essentially because 
non-linear interaction terms in the Hamiltonian easily make the problem chaotic and 
intractable. Only very few models, such as the harmonic oscillator, are solvable exactly. 
Somewhere in between these two extremes live the integrable models. They may contain 
highly non-trivial interactions, but they can nevertheless be solved completely by the 
right ansatz for the wave functions. Such a wave function will depend on a couple 
of parameters, and quantization conditions will impose a system of equations on them. 
There need not be a general analytic solution to these equations — after all the spectrum 
of integrable models is usually highly non-trivial — nevertheless the reduction to a small 
number of parameters is sufficient to make the spectral problem much more tractable 
than for generic quantum mechanical models. 

Typical integrable models are formulated in 1+1 or 2 dimensions. They include field 
theories, non-linear sigma models, particle models, vertex models and spin chains (we 
shall focus on the latter in this article). A plethora of integrable models with all kinds of 
features is known to date, and it appears near impossible to make a complete census. A 
central insight towards this goal was made in the 1980's by the Leningrad/St. Petersburg- 
School who related integrability to the existence of large hidden symmetries. Through 
the enumeration of suitable symmetry algebras one can hope to classify the integrable 
models. 

For example, a very large class of integrable spin chains can be derived and in- 
vestigated using Yangians and quantum affine algebras U g (g) |1|. In particular, the 
Heisenberg XXX (algebraic) spin chain and its relatives with different symmetry algebra 
and/or different representation are all based on the Yangian double of q. Likewise, 
the quantum-deformed XXZ-like (trigonometric) spin chain cousins are related to the 
quantum affine algebra XJ q (g). XYZ-like (elliptic) spin chains also have a similar but 
much more elaborate underlying symmetry algebra. 

A famous integrable spin chain model that has escaped this classification for a long 
time is the one- dimensional Hubbard model, see |2|. It is a model of electrons propagating 
on a chain of nuclei. Each nucleus site can either be unoccupied, singly occupied with 
electron spin pointing up/down or doubly occupied with opposing electron spins. In 
total there are four states per site (c' a is a fermionic electron creation operator) 

1°), IT)~4|o>, u>~4|o>, it)~cl4|o). (i.i) 

The middle two states are considered fermionic while the outer two states are overall 
bosonic. This model is exciting because it shows some characteristics of superconduc- 
tivity, and therefore it is very desirable to understand its foundations well. Integrability 
was established by Lieb and Wu who also solved the spectrum by means of the Bethe 
ansatz 3 1 . An R-matrix which encodes the integrable structure was later found by Shas- 
try |4|. On the one hand, the R-matrix is the foundation for much of the integrable 
machinery, such as the algebraic Bethe ansatz |5,6|. On the other hand, this particular 
R-matrix is rather exceptional because unlike most other known R-matrices it cannot 
be written as a function of the difference of two spectral parameters. Altogether, the 
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algebraic origin of the R- matrix remained mysterious. It is well-known that it is symmet- 
ric under two undeformed su(2) algebras: spin and (twisted) eta-pairing |7| symmetry. 
Therefore one may expect the underlying algebra to be of Yangian (algebraic) rather 
than of quantum affine (trigonometric) type. Indeed, two su(2) Yangians algebras were 
identified |8|, but they are not sufficient to explain the R- matrix. A fusion procedure of 
two XX models was used to derive the R-matrix and explain its features, but it seems 
very specialized to the model at hand and it hardly illuminates the symmetries. 

New insight into the algebraic structure came from a very different and unexpected 
direction: the field of gauge theory, string theory and the so-called AdS/CFT correspon- 
dence, which relates certain pairs of gauge and string theories. In that context it was 
observed that M = 4 superconformal Yang-Mills theory in the 't Hooft limit and its dual, 
IIB string theory on the AdS^xS 5 , both display signs of integrability |9|, see |10.11| for re- 
views and further references. The asymptotic coordinate Bethe ansatz |12| for the gauge 
theory spin chain leads to (two copies of) an interesting scattering matrix |13| which is 
not of difference form. A construction for strings in light-cone gauge |14| leads to an 
equivalent S-matrix |15|. The scattering particles have four flavors {(ft 1 ) , \(fi 2 ) , {ip 1 ) , \tp 2 ) , 
the former two being bosonic and the latter two being fermionic. The set of particle 
flavors is equivalent to the states of a site in the Hubbard model 

|o)~|^), |T>H^>, II) ~hA 2 >, IDH0 2 >, (i-2) 

and it was observed that the S-matrix has a structure reminiscent of the R-matrix for 
the Hubbard model |16|. Indeed, the two matrices can be mapped into each other 
exactly |17| which leads to a very curious connection between string theory and the 
integrable structure of the Hubbard model. This link is also reflected in the asymptotic 
Bethe equations for planar AdS/CFT |18| which contain (two copies) of the Lieb-Wu 
equations in disguise. 

The large amount of supersymmetry present in the string/gauge theory system, the 
superalgebra psu(2,2|4), thus made its way into the integrable structure of the one- 
dimensional Hubbard model: What remains of this symmetry in the above scattering 
picture is (two copies of) su(2|2) |10|. A crucial point for the further understanding was 
that the symmetry is centrally extended by gauge transformations inherent to the gauge 
theory |13| or by residual transformations in the light cone gauge for string theory |19|. 
The symmetry of the S-matrix turns out to be an (exceptional) threefold central extension 
P) of the psu(2|2) superalgebra 1 

P) := psu(2|2) x M 3 = su(2|2) x R 2 . (1.3) 

This algebra contains the two well-known bosonic su(2) symmetries of the Hubbard 
model which relate the two bosonic and two fermionic states, respectively. However, 
the additional fermionic generators of the algebra also relate the bosons to the fermions 
and vice versa. In fact the algebra is strong enough to fully constrain the form of the 
R-matrix |13|. 

1 For simplicity we shall consider the algebra to be complex and do not distinguish between psu(2|2), 
psl(2|2) or pst(2|2, C). Reality conditions refer to the real version psu(2|2) of the algebra. 
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The proper framework for the symmetries of the R-matrix and thus for the integrable 
structure of the one- dimensional Hubbard model is expected to be a quasi-triangular 
Hopf algebra |20|. The goal is then to find the universal R-matrix of which the above 
R-matrix is the fundamental representation. However this requires to first identify the 
complete symmetry algebra of the R-matrix. Generically one may expect the algebra 
to be a Yangian double: a deformation of the universal enveloping algebra of the loop 
algebra of the underlying symmetry f). Indeed many of the Yangian generators have been 
identified |21,22|. An investigation of the classical limit |23| of the R-matrix has then 
revealed the complete classical structure in terms of a quasi-triangular bialgebra |24 
based on a curious deformation of the loop algebra u(2|2)[u, u" 1 ). This result shows that 
in addition to the central charges, there must also be inner automorphisms. 

At least two important steps remain to be taken: First, the classical bialgebra needs 
to be quantized to a Yangian double. Second, the universal R-matrix for the Yangian 
needs to be established which makes the Yangian into a quasi-triangular Hopf algebra. 
However, it is not easy to deal with Yangian doubles and their algebraic structure because 
proper quantization of the higher levels is somewhat unintuitive and specialized to the 
algebra g. Instead one usually considers the corresponding quantum affine algebra U g ($j) 
of which the Yangian is a contraction limit for q — > 1. Here one pays the price that the 
q symmetry is quantum deformed and not as manifest as in the Yangian. Instead one 
gains a uniform treatment for the quantum deformation of the Kac-Moody structure of 
the whole of the affine algebra q. It is the aim of the present paper to lay the foundations 
for the quantum deformation of the integrable structure of the one-dimensional Hubbard 
model. Here we shall start with the quantum deformation U 9 (f)) of the algebra f) and 
leave the full quantum affine algebra U g (fj) (or rather the deformed U g (u(2|2)[w, u^ 1 ])) 
and its quasi-triangular structure for future work; the corresponding Yangian double 
would follow as the limit q — > 1. We will then derive the fundamental R-matrix which 
should be understood as the quantum deformation of Shastry's R-matrix for the Hubbard 
model. We apply it to derive the Bethe equations for periodic wave functions and a three- 
parameter family of Hamiltonians with U g (su(2) x su(2)) symmetry which includes the 
Hubbard Hamiltonian as a special case. 

In fact, many attempts have been made to modify and generalize the Hubbard Hamil- 
tonian due to the exceptional properties of the Hubbard model. Widely discussed mod- 
ifications are the EKS model |25|, the supersymmetric U- model |26| and its quantum 
deformation |27| as well as the su(n) Hubbard models |28| These can all be explained 
with the available integrability toolkit: The EKS model is a model based on the fun- 
damental representation of u(2|2) and the supersymmetric U-model is based on the 
four-dimensional representation of su(2|l). These models are somewhat similar to the 
Hubbard model, but they do not include it as a special case. The su(n) Hubbard models 
employ an external coloring of states which preserves integrability. This coloring can be 
applied to any integrable model with conserved charges |29| and it does not alter the 
underlying symmetry. Further similar models have been discussed in |30|. An important 
class of models which is also discussed in this context consists of the supersymmetric t-J 
model 1 31 1 and some of its deformations such as the Bariev model |32| and others |33|. 
The main difference is that these models use a three-dimensional representation on each 
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site and thus the Hilbert space is very different from one of the Hubbard model. The only 
known true deformation of the Hubbard model appears to be a Hamiltonian composed 
by Alcaraz and Bariev |34|. The Hamiltonian contains substantially more terms and so 
far it has not been investigated further in the literature. The Bethe equations for this 
model were given in |34|, and the bear some resemblance with those for the XXZ model. 
Thus it is conceivable that the Alcaraz-Bariev model is a quantum-deformation of the 
Hubbard model. We shall address the question whether we can recover this Hamiltonian 
at the end of our work. 

The present paper is organized as follows: We start with a technical part concerning 
the algebra, R-matrix and Bethe ansatz in Sec. 2 3 4, respectively. In the second part 
consisting of Sec. 5 we apply the obtained results to a concrete spin chain model. It is 
not necessary to read the earlier sections (in full detail) to understand the later sections. 

First the Uq(f)) symmetry is introduced in Sec. 2. We also sketch finite representations 
of the algebra with particular focus on the fundamental representation needed for the 
derivation of the fundamental R-matrix in Sec. 3, We then perform the nested Bethe 
ansatz for this R-matrix in Sec. 4 to obtain the Bethe equations for a periodic chain. 
Finally, in Sec. 5 we derive a class of integrable Hamiltonians associated to the R-matrix. 
These constitute quantum deformations of the Hubbard Hamiltonian. In particular, we 
recover one of the models proposed by Alcaraz and Bariev. We conclude and give an 
outlook of open problems in Sec. 6. 

2 The Hopf Algebra U q (su(2\2) x R 2 ) 

We start with the the construction of the symmetry algebra underlying quantum defor- 
mations of the one-dimensional Hubbard model. This part is rather technical in nature 
and can be skipped or be used as a reference for the following sections. 

2.1 From su(2|2) to U(su(2|2) k R 2 ) 

We start by introducing the universal enveloping algebra U(su(2|2) x IR 2 ) step by step 
starting from the Lie superalgebra su(2|2). This will help us to understand and derive 
the quantum deformation U g (su(2|2) k M 2 ). 

Lie Superalgebra. The Lie superalgebra su(2|2) is generated by the su(2) x su(2) 
generators SH a b, £ a /3, the supercharges a &, & a /3 and the central charge £r The Lie 
brackets of the su(2) generators take the standard form 

[<R\ e c 5 ] = 5 c b e a 5 - \5 a b & 5 , G c s ] 

2 To obtain the simple Lie algebra psu(2|2) from su(2|2) we would have to project out this central 
clement. 



= -5^ p + \5^ 5 . (2.1) 
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Figure 1: Distinguished Dynkin diagram of su(2|2). 



The Lie brackets of two supercharges yield 

{in, & c 5 } = 5 c b £ a 5 + 5%ft c b + 5 c h 5tt. (2.2) 
The remaining Lie brackets vanish. 



Central Extension. This algebra has two further possible central extensions ^3, .ft. 
They are generated by Lie brackets of alike supercharges 

{£A, £T d } = e^e bd ^, {6V & c s } = e ac e p5 Si. (2.3) 

The centrally extended algebra with these charges shall be denoted by 

h :=5u(2|2) xM 2 = psu(2|2) xM 3 . (2.4) 



Universal Enveloping Algebra. The universal enveloping algebra U(h) of f) is gen- 
erated by polynomials of the Lie algebra generators. The Lie brackets are represented 
as commutators or anti-commutators (depending on the statistics of generators) 

[X, Y] -> XY - YX, {X, Y} — > XY + YX. (2.5) 

They respect the Lie algebra relations by identification of certain polynomials, e.g. 

a 6 £F d + 0.\£l a b = e ai e bd ^. (2.6) 



Chevalley Basis. Within the universal enveloping algebra it is not necessary to keep 
all generators of the Lie algebra explicitly. For example, the central charge ^3 can be 
represented through a quadratic combination of supercharges, see (2.6). A minimal set 
of generators for this rank-three algebra is given by three Cartan generators Sjj, three 
simple positive roots <Bj and three simple negative roots j = 1, 2, 3. One may identify 
them with the Lie generators as follows 





= m 2 2 - 


m\ 


= 2<K 2 2 , 


Ci 


= ^ 2 i, 




= JH 1 2 




= -€- 








= n 2 2 , 




= 6 2 2 


Sd3 


= £ 2 2 - 




= 2£ 2 2 , 


<S 3 


= #2, 




= £ 2 i- 



This basis corresponds to the distinguished Dynkin diagram of bu(2|2) in Fig. 1. The 
symmetric Cartan matrix in this basis reads 

/ +2 -1 \ 
A jk = -1 +1 . (2.8) 
V +1 -2 / 
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Note that the algebra su(2|2) has a degenerate Cartan matrix and therefore has a null 
vector which reads Vj = (1, 2, 1). It can be used to express the central charge of su(2|2) 

1 3 

€ = "a £ v & = - " «2 " 1*3- (2.9) 
i=i 

The other two central charges in the Chevalley basis take the form 

¥ = { [<*! , <E 2 ] , [<£ 3 , <£ 2 ] } , £ = { [Si , S2] , [5s, S2] } • (2.10) 

Commutation Relations. Let us now state the commutation relations of the uni- 
versal enveloping algebra. Commutators with the Cartan generators S)j are given by 
(j, k = 1,2, 3) 

[%£ fe ] = 0, [% <B k ] = +A jk <B k , [SjM = -A jk $ k . (2.11) 
The non-trivial commutators of positive and negative simple roots read 

[<£i,Si]=+Si, {£2,^2} = -*2, [€3,^3] = -«3. (2.12) 

Note that we have normalized the generators in a way such that all relations can be 
expressed using the symmetric Cartan matrix Aj k . For superalgebras this leads to a neg- 
ative sign in [<£j, Sj] for one of the two bosonic subalgebras. The remaining commutators 
between positive and negative simple roots vanish in this basis 

= ioij^k. (2.13) 

Finally, we need to impose the Serre relations between positive and between negative 
simple roots. 

= [(*!, £ 3 ] = £ 2 £ 2 = [<£i, [<Si, £2]] = [(£ 3 , [€ 3 , <£ 2 ]] 
= [Si, 5a] = S2S2 = [Si, S2]] = [S3, [S3, S2]] • (2.14) 

Note that for superalgebras these standard Serre relations are not sufficient. For the 
algebra U(su(2|2)) we need two additional Serre relations which effectively read 

= = (2.15) 

It is however consistent to drop them altogether which leads to the centrally extended 
algebra U(f)). 



2.2 Quantum Deformation 

The quantum algebra U q {g) is a deformation of the universal enveloping algebra U(g) 
of some Lie algebra g. It is obtained by inserting various factors and exponents of q 
in various relations. A convenient combination to use in this context is the so-called 
quantum number 

[ X]q :=^£l. (2.16) 
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This relation is defined in the same way if x is a generator. There are two ways to achieve 
this: Either one defines q x = 1 + x log q + \x 2 log 2 q + . . . as a formal power series. For 
practical purposes one would assume that q « 1 and thus log q ~ 0. Alternatively one 
can define q x as an element of the Hopf algebra and q~ x as its inverse. This is particularly 
useful if x appears only with integer multiplicity in exponents, q nx = {q x ) n , as will be 
the case here. 

In the previous section we have seen that the two central charges can be obtained by 
dropping two Serre relations. It is straightforward to apply the same central extension 
for the quantum deformed algebra U 9 (su(2|2)) to obtain U g (f)). 

Commutation Relations. Let us begin with the deformation of the commutation 
relations of the algebra U g (su(2|2)). The commutators with Cartan elements do not 
receive deformations 

[%#fc] = 0, [fij, <£ k ] = +A jk £ k , [%&] = -AjkSk. (2.17) 

The Cartan elements usually appear in exponents, and it is convenient to note the 
exponentiated form of these relations 

q*lE k = q +A *<£ k q^, q%$ k = q~ A ^ k q^ , (2.18) 

The non-trivial commutators of simple roots read in the deformed algebra 

[<Si,3l] = [f>i] q , {£ 2 ,S 2 } = -[f>2] q , [<S 3 ,&] = -[3 3 ] 9 , (2-19) 

and the remaining mixed commutators vanish 

= hij^k. (2.20) 

The Serre relations have the same form as in the undeformed algebra (2.14), but with 
some additional factors of q due to the quantum- adjoint action. Spelled out they yield 

= [(£!, <S 3 ] = [ffl, &] = «2<£ 2 = fofo (2.21) 

= £i£i<£ 2 -(q + q- 1 )^^! + € 2 Ci€i = <S 3 <B 3 <£ 2 ~(q + g" 1 )^^ + £2 £3^3 

= $1$1$2 - (q + q~ X )$vSl!S\ + $2$X$\ = fofo&Z ~{q+ <7~ 1 )#332&5 + #2#3#3- 

Central Elements. The standard central element from the Cartan subalgebra remains 
undeformed 

£ = - $2 - §$3. (2.22) 

As before, the ordinary Serre relations obtained from the Cartan matrix are not suffi- 
cient to define U 9 (su(2|2)), but we need the constraint ^3 = ^ = with the quantum 
deformation of (2.10) 

<p = £i£ 2 £ 3 £ 2 + £ 2 £ 3 £ 2 <£i + ^2^2 + £ 2 £i<£ 2 £ 3 - (5 + 5 -1 )€ 2 €iC3^2, 

£ = Sifo&sfo + fofofoSi + fofofofo + fofofofo - (<? + g~ 1 )3 r 2T$i3 r 33 r 2- (2-23) 

In the centrally extended algebra U g (h) we will not impose the constraint and thus obtain 
two non-trivial central elements ^3, It is straightforward, if tedious, to confirm that 
these polynomials are indeed central elements of the quantum-deformed algebra. 
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2.3 Representations 

It is commonly the case for a finite-dimensional simple Lie algebra g and for generic values 
of q that the representation theory of the quantum deformed algebra U g (g) is analogous 
to the one of q. 3 We have studied some of the simplest finite-dimensional representations 
of Uq(f)) and found agreement with this rule. Here we would like to give an overview 
over some basic representations of U g (f)) assuming that there is no qualitative difference 
to the undeformed case. 

Outer Automorphism. Finite representations of f} were studied in |17|. The starting 
point was the representation theory of su(2|2) with fixed eigenvalue Co of the central 
charge £, studied in, e.g. |35,36|. The SL(2) outer automorphism of 1) then allows to 
rotate the triplet of eigenvalues (C , 0, 0) of the central charges (£, ^3, &) to any desired 
triplet (C,P,K) with C5 2 = C 2 — PK. 4 The representation of the psu(2|2) generators 
0,(5 is obtained by conjugation with the SL(2) element. Note that the combination 
£ 2 = £ 2 — is invariant under the automorphism. 

It appears that for the quantum-deformed algebra U q (t)) there also exists a simi- 
lar outer automorphism. We could use it to relate representations of U g (f)) to those 
of U g (su(2|2)), see |37|, which in turn are analogous to those of su(2|2). We how- 
ever do not yet understand the automorphism explicitly and therefore the existence of 
the below representations is an educated guess. The above combination €? should be 
quantum-deformed to some [£ 2 ] g . The classification of representations would then use 
the eigenvalues of this operator. We find that the smallest representations are indeed 
characterized by the eigenvalues of the operator 

[l\:=[<Lf q -^. (2.24) 

Presumably this combination is invariant under the tentative automorphism. 

Let us consider typical (long) and atypical (short) representations of U g (f)) in analogy 
to the su(2|2) representations studied in |35|. 

Long Multiplets. The standard finite-dimensional type of representation shall be 
denoted by 

{m,n;C,P,K}. (2.25) 

It corresponds to the typical highest-weight representations of su(2|2) with Dynkin labels 
[m; r; n) and r = ±([C 2 ],j) 1//2 — \n + ^m. The non-negative integers n,m represent the 
Dynkin labels of the su(2) x su(2) subalgebra. 

We can decompose a U g (P|) multiplet into irreducible multiplets of the subalgebra 
Uq(su(2) x su(2)). Let the symbol [k] represent the U g (su(2)) representation with spin 

3 In quantum algebras one singles out the case when q is a root of unity because the representation 
theory is very special at these points. 

4 For simplicity we shall consider the algebra \) to be complex, and thus there is no distinction between 
positive and negative values of C 2 — PK. 
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fc/2; [—1] is the empty (zero-dimensional) representation. Then the long multiplet de- 
composes as follows 

{m,n} -> ([m] © [1] © [1], [n]) © ([m], [n] © [1] © [1]) 

© ([m] © [1], [n] © [1]) © ([m] © [1], [n] © [1]). (2.26) 

Note the well-known tensor product of U 9 (su(2)) representations 

min(m,n) 

[m] © [n] = [m + n-2Jfe]. (2.27) 

fc=0 

The dimension of a long multiplet is thus 16(m + l)(n + 1). 

Short Multiplets. A short multiplet shall be labelled by 

(m,n;C,P,K). (2.28) 

It can only exist when the constraint 

[C\ = [Cf q - PK = [\{m + n + l)}l (2.29) 

holds. The decomposition into irreducible representations of U g (su(2) x su(2)) takes the 
form 

(m,n) -> ([m] © [1], [n]) © ([m - 1], [n - 1] © [1]) 

© ([m], [n] © [1]) © ([m - 1] © [1], [n - 1]). (2.30) 

This multiplet has dimension 4(m + l)(n+ 1) + 4mn. For example, the four-dimensional 
fundamental representation is the special case m = n = which decomposes into (0, 0) — ► 

([i],[o])e([o],[i]). 

Multiplet Splitting. The vector space of a long multiplet {m, n; C, P, if} can be 
decomposed into the vector spaces of two short multiplets 

{m, n; C} -> (m + 1, n; C 5 ) © (m, n + 1; d). (2.31) 

However, the typical representation of U g (fj) is in general irreducible. A special case is 
when the labels and central charges obey the shortening condition 

[C\ = %{m + n + 2)]% (2.32) 

which is equivalent to the constraint (2.29) for the short multiplets in (2.31). The long 
representation may then be reduced into the above two short representations. However, 
in general one cannot expect the long representation to be decomposable, it merely has 
one short subrepresentation which closes on one of the smaller vector spaces. When 
projecting out this small vector space one obtains a short factor representation on the 
other small vector space. This fact is perhaps best illustrated by a figure, see Fig. 2, 
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Figure 2: Multiplet Splitting. The long multiplet {m,n} consists of two short multiplets 
(m + l,n) and (m,n + 1) and a representation of U g (f)) maps between the short multiplets 
(arrows). Generically, the short multiplets are connected in all possible ways (left). When the 
shortening condition holds one (middle) or both (right) arrows between the short multiplets 
are broken. In the second figure from the left, (m + 1, n) is a subrepresentation while (m, n + 1) 
is a factor representation, and the long multiplet is indecomposable. The long multiplet in the 
right figure is fully decomposable. 

Tensor Products. Let us state some formulas for the decomposition of tensor products 
of long and short representations. It is convenient to introduce some short hand notation 
for sums of similar representations 

{[k] © [l],n} := {k,n} © {l,n}, etc.. (2.33) 

The tensor product of two long representations of U 9 (f)) decomposes as follows 

{m,n} © {k, 1} = {[m] © [k] © [1] © [1], [n] © [/]} 
© {[m] © [k] © [1], [n] © [Z] © [1]} 
© {[m] © [k] © [1], [n] © [/] © [1]} 
© {[m] © [k], [n] © [Z] © [1] © [1]}. (2.34) 

For a tensor product of a long and a short representation the following decomposition 
applies 

{m,n} © (k,l) = {[m] © [k] © [1], [n] © [/]} 

© {[m] © [fc - 1] © [1], [n] © [Z - 
© {[m] © [A;], [n] © [Z] © [1]} 
© {[m] © [A; - 1], [n] © [Z - 1] © 

Finally, the tensor product of two short representations reads 

(m,n) © (k,l) = {[m] © [k], [n] © [I]} 

©{[m]©[A;-l],[n]©[/-l]} 
©{[m-1]® [A;],[7i-1]® [Z]} 
© {[m - 1] © [fc - 1], [n - 1] © [/ - 1]}. (2.36) 

These formulas agree with the decompositions (2.26 2.30). 
2.4 Hopf Algebra 

Here we will complete the Hopf algebra structure of our quantum deformed algebra. We 
first define the coproduct and then state the unit, counit and antipode. 



1]} 

[1]}. (2.35) 
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Coproduct. We shall use the standard coproduct for the generators of U 9 (su(2|2)) 

A(l) = 1(8)1, 
A(Sjj) =$jj <8> 1 + 1®% 
A((Bj) = (Bj <g> 1 + <T^ <g> Cj, 

A(^)=3 r J -®g fl, ' + i®3 r j, (2.37) 

where 1 denotes the unit element of the Hopf algebra. Note that the relations are fully 
compatible with the Serre relations and commutators. The coproducts for the central 
charges follow by substituting their expressions in terms of the above generators 

A(C) = £® 1 + 1 (8)£, 
A(«p) =<p®l + g 2£ ®<p, 

A(£) = £® g" 2C + 1 <8>£. (2.38) 

Note that the coproduct for the generators ^3, ^ turns out to be proportional to the 
generators themselves as expected. This fact is related to the the role of = A = as 
Serre relations of U ? (su(2|2)). The coproduct of such consistent identifications has to be 
proportional to these or other identifications. 

Unit, Counit and Antipode. The unit element of the Hopf algebra is denoted by 
77(1) = 1. The counit e : U g (f)) — > C takes the usual form 

e(l) = 1, eft;) = e(fBi) = e(&) = 0. (2.39) 

Finally, the antipode S : U g (f)) — > U g (f)) is uniquely fixed by the compatibility condition 

H o (S <g> 1) o A(3) = a* o (1 (8) 5) o A(3) = ?7 o e(3), (2.40) 

for all 3 e U q {i)), where 77 : C — > U g (h) denotes the unit and fi : U 3 (f)) (8> U g (f)) — > U 9 (f)) 
denotes the product. One gets 

5(1) = 1, S(^) = -% 5(6,-) = = (2-41) 

and similarly for the central charges 

5(e) = 5(<p) = -<T 2C< P, S(£) = -q 2€ &. (2.42) 

Hermitian Conjugation. So far we have for simplicity assumed a complex algebra. 
To restrict to a real algebra we must identify the generators with their conjugates. The 
proper hermiticity relations compatible with the coproduct are 

f>) = % e] = q %v,. (2.43) 

This implies the following relations for the central charges 

£t = £ ? qjt = ? 2££ 

To have a hermitian coproduct we furthermore have to constrain q to be real. 
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Braiding. Just as for the undeformed algebra U(f)) we shall deform the coproduct 
slightly and according to a Z-grading of the algebra. This grading associates the charges 
+2, +1, —1, —2 to the generators ty, (£ 2 , #2, respectively; the other generators are un- 
charged. This braiding will later lead to a very non-trivial R-matrix. 

For the braiding we introduce a new abelian generator il and deform the coproduct 
(2.37) as follows 

A(<5 2 ) = <£ 2 <8> 1 + g-^il <8> £ 2 , 

A(ff 2 )=ff 2 ®g fla +a- 1 ®^ 2 , 
A(«P) = <P® l + g 2£ il 2 ®<p, 
A(JS) = £ g> <T 2e; + iT 2 g> 

A(il)=il®il. (2.45) 

The coproduct of the other elements remains undeformed. This deformation of the 
coproduct is consistent with the algebra because the exponents of il follow the Z-grading 
of the algebra. The counit for il should be defined as e(il) = 1. The antipode must then 
obey the relation (2.40) which leads to the following modifications of (2.41) 

S(<E 2 ) = -f^iT 1 ^, S(ff 2 ) = -<T fe il£ 2 , 5(il)=il- 1 (2.46) 

and for the central charges 

S(£) = S($S) = -q- 2€ il- 2 ¥, S(R) = -q 2€ ii 2 &. (2.47) 

The hermitian conjugation remains untouched by the braiding and 

il f = il -1 . (2.48) 

Cocommutativity of the Center. A Hopf algebra is called quasi-cocommutative 
if the coproduct A (3) and the opposite coproduct A op (3) := V(A($)) are related by 
conjugation. Here V defines the graded permutation operator for the tensor product. The 
conjugation is specified by an R-matrix 1Z G U g (f))®U g (f)) satisfying the cocommutativity 
condition 

A op (Z)1l = nA(Z). (2.49) 

For generators 3 from the center of U g (fj) this relation implies that the coproduct 
must cocommute A op (3) = A (3). With our coproduct this is obvious for the central 
charge £ 

A(£) = £®1 + 1®£ = Aop(£). (2.50) 
The other two central charges ^3 and & however do not enjoy this property 

AC£)=<p®l + g 2£ il 2 ®<p, A op (^)=^®ilV c + l®^, 

A(£) = 8.® q- 2£ + ii~ 2 (g)&, U A op (£) = &(g)!d- 2 + q- 2€ ®&. ^' ' 

By taking the difference of the coproducts one can see that the central charges ^3, & have 
to be identified with € and il as follows 

<P = ^ a (l-g 2C il 2 ), & = gaT 1 (q- 2C -SX- 2 ). (2.52) 
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Here g and a are two global constants of the reduced algebra. 

An alternative derivation of the braiding and identifications based on a consistent 
factorised scattering picture is presented in App. B. 

2.5 Fundamental Representation 

The fundamental representation is the short multiplet 

(C,P,K) = (0,0;C,P,K). (2.53) 

The vector space is spanned by four states 

l^ 1 ),^ 2 ),^ 1 ),^ 2 ), (2.54) 

where \<p a ) are bosonic and \ifj a ) are fermionic. The fundamental action of the Chevalley- 
Serre generators is given by 

iD 1 |0 1 ) = -|0 1 ) ) = -{C-±W), e 1 |0 1 ) = g+ 1 /2| 2 ); £ 2 | i ) = f# i }; 

^|0 2 ) = +|0 2 ), Sy 2 \f) = -(C+\)\<p), £ 2 |0 2 > = a|^ 2 >, &|^> = ?- l/ V>, 

w 2 ) = +\r), s)2m = -{c+\w), £ 3 m = q - i/2 m, &iv 2 > = <# 2 >, 

(2.55) 

The braiding generator il will act with the same eigenvalue U on all four states. 

Constraints. It is not too hard to see that the closure of the algebra of supercharges 
leads to the constraints 

ad=[C+l] q , &c=[C-§]„ ab = P, cd = K. (2.56) 

This is in agreement with the constraint (2.29) which is required for the fundamental 
representation 

[C\ = W\l (2-57) 
It furthermore follows that the constraint ad — be = 1 for q = 1 is deformed to 

{ad - qbc)(ad - q~ l bc) = 1. (2.58) 

Note that there is a subtlety in taking the limit q — > 1 in the above constraint (2.58): 
It leads to two inequivalent constraints ad — be = ±1 which correspond to different 
embeddings of the Lie algebra into the Chevalley basis. The plus sign yields the above 
constraint with the identifications in (2.7) or equivalently (E 2 = — H 2 2, ^2 = — © 2 2- The 
minus sign corresponds to one of the following embeddings 

£ 2 = tH 2 2 , S 2 = ±© 2 2 , £2 = t + ifr + AS 3 - (2.59) 
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x^ 1 Parameters. For later convenience we shall introduce the new representation pa- 
rameters x , 7 and rewrite a, b, c, d as follows 

a = y/gnf, 

b=y^ — (x--q 2C - 1 x + ), 

7 x~ 



iy/gi q 



-C+l/2 



a 

d= hH q c+i/2 f x - _ q -2C-i x+ ^ {2m) 

In terms of these parameters the constraint (2.56) implies the following quadratic relation 
between x ± 

x + q 1 . i . / x + qx~ \ i , 

— + ^L-qx- +zg ( q - q -l)l =_. 2.61 

q x + qx~ \qx~ x + J g 

The central charge C cannot be written unambiguously using x ± , but the combination 
q 2C is well-defined 

n 2c _ „ (g - q~ l )/x + - ig' 1 _ x (g - g- 1 )^ + ig- 1 
\q~q )/ x — % 9 \q — q ) x + % g 

These two expressions are equivalent upon (2.61). Finally, the central charges P, K read 

P = ga(l- q 2C ^ , K = ga- 1 [q- 2C - . (2.63) 

With the above identification (2.52) of and .ft with il, the squared eigenvalue of the 
latter on the fundamental representation reads 

U 2 = — . (2.64) 
qx 

Fermion Normalization. The parameter 7 adjusts the normalization of fermions \ip a ) 
with respect to bosons \<f) a ). Furthermore the parameter a adjusts the normalization of 
(£ 2 , vs. #2, A particularly useful choice for 7 will turn out to be 

= y/-iaq c+1 / 2 U(x+ -x-) 

Apart from the features discussed below, we expect it to have nice analytic properties 
analogous to the case of the undeformed case discussed in |38|. Despite some simplifica- 
tions we will largely keep a and 7 unspecified in this paper. 
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Uniformization. For a fixed parameter g the constraint (2.61) defines a complex torus 
with modulus 

k = teg^Jl - g 2 {q - q' 1 ) 2 . (2.66) 

The quotient of the parameters x can be expressed conveniently through the point z 
on a doubly-periodic complex plane using 

x + cn(z) +irsn(z) 

(2.67) 



qx~ cn{z) — ir sn{z) 
where the constant r is given by 



/ ! + 16^-16^-^-1)2 
( 1/ l+V(gV2 + g -i/2)2 • 

The individual parameters x ± are then written as rational functions of sn(z), cn(^), dn(z), 
but we shall refrain from stating these functions explicitly here. 

Transposition. The antipode map corresponds to transposing the representation ma- 
trices 

S(3) ~ C~ X TC, (2.69) 

where the supertransposition is defined by A S A. := (—l)^ +1 ^Akj and C is the charge 
conjugation matrix 

C\<f) = -g +1 / 2 |0 2 ), C|^> = -g +1 /V>, 

C|0 2 ) = +g^/2| i ); C |^2 ) = +9 -i/2|^,i). ^ > 

The conjugation matrix was constructed such that the equation is satisfied for the bosonic 
Chevalley-Serre generators (Si, (£3, #1, $3. By solving this equation for the fermionic 
generators £2, $2, we obtain the antipode representation parameters a, b, c, d in terms of 
the original parameters a, b, c, d (see also |15|) 

a = -U-'q-^b, 
b = U- 1 q- c - 1 '\, 
c = -Uq c+1 ' 2 d, 

d = Uq c - 1/2 c. (2.71) 
We can solve these equations for the crossed spectral parameters and get 5 

a 1/2 a 

x+ = s(x + ), x- = s(x~), 7 = - (q c U -q-Clf- 1 ). (2.72) 

7 

Note that 7 is compatible with the choice (2.65) of 7 as a function of x . The antipode 
map s(x) for the parameters x is defined as 

ig 1 x — [q — q 1 ) 



5 Note that 7 = —7 which is consistent with the Z/j-periodicity of supertransposition. Alternatively 
one could define U = —l/U leading to 7 = 7. 
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This map x i— > s{x ± ) acts on the braiding factor and the central charge as follows 

C=-C, U=±. (2.74) 

Hermiticity. Let us introduce the canonical scalar product for our states 

(<P a \<P b ) = 5 b a , (Vg^>=e (2.75) 



The representation is hermitian if the conditions in (2.43) hold. They imply the relations 
a* = q c+l l 2 d and b* = q c ~ 1 ^ 2 c. Expressed in terms of x ± -parameters they correspond to 

(0* = -At- (2-76) 



s(x~) s(x~ 



This is equivalent to the condition that the uniformized parameter z on the complex 
torus is real. Furthermore the moduli of the parameters 7 and a are constrained 



| 7 | 2 = -i x + +iq 2 °+ 1 x - , \a\ 2 = l. (2.77) 

Again, the first constraint is automatically satisfied when 7 is given by (2.65) as a function 
of x . 



3 The Fundamental R-Matrix 

To determine the full universal R-matrix is a formidable task. Here we will satisfy 
ourselves with the lesser task of finding the representation of the R-matrix on two fun- 
damentals. 



3.1 Matrix Structure 

We construct the fundamental R-matrix by demanding that it satisfies the cocommuta- 
tivity relation (2.49) 

A op (z)n = nA(z). (3.1) 

The most general ansatz which satisfies the relation for the U g (su(2) x su(2)) subalgebra, 
i.e. for £1, £3, #1, #3, is given in Tab. 1. The ten functions A, B,C, D, E, F,G, H, K,L 
are a priori free. They are however fully constrained up to one overall phase factor R\ 2 
by cocommutativity w.r.t. the supersymmetry generators (£2, #2- This finding is related 
to the irreducibility of the tensor product of two fundamental representations (2.36) 

(0,0) ® (0,0) = {0,0}. (3.2) 

We present our findings for the ten functions in Tab. 2, This derivation parallels com- 
pletely the case of the undeformed algebra in |17| and for q = 1 the results agree. Note 
that C12 = F12 and H 12 = K 12 if 7 is determined by (2.65). 
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n\ 
n\ 



a ii 



,li2\ 



n\<j> 2 <j> 2 ) 



12 



£.1/1 



qA 12 + q 1 B 12 



(.2/1 



g + g -1 



.4 



12 



g + g 



-^ 12 lJ.lj.2\ ? 1 Cl2| / 2,l\ , ^12 I ,1 /2\ 

g + g 1 g + g 1 



.4 



12 



B 



^y>+ * ?":_r i! w > 



g M12 + g-Bi 2| ^ 



g + g- 



C 



12 



T IV>V> 



g + g 1 



g + g- 



gC 



12 



g + g- 



A 12 |0 2 2 ) 
,i„/,i 



£>1 2 |V>V> 

gL>i2 + q~ x E 12 



g + g 1 
D 12 -E 



D 



12 



"^12 I / 1 1 2\ , 9 ^ 



12 | ,2,1 



g + g 

-A2|^V) 



g + g" 

— 12 I / 2 / 1\ « ^12 -T y-^12| ,l ,2\ 

— |v> V ) — „ , „_ 1 — f V» ) 



g + g£i2 



g + g 1 
F12 



12 



g + g 



g + g 



=rlW>4 



g + g _1 

qFi 2 
q + g _1 



G 12 |0°^} + ^i2|^0 a ) 
i^V^ + WV} 



-U2\ 



Ui2\ 



Table 1: The fundamental R-matrix of U 9 (fj). 



3.2 Discrete Symmetries of the R-matrix 

Braiding Unitarity. An R-matrix is expected to obey the so-called unitarity condition 

TZ 12 TZ 21 = 1 ® 1. (3.3) 

This condition says that the operation of flipping the order of two sites along a chain 
is an involution. Unitarity implies the following ten relations for the coefficients of the 
operator in Tab. 1 

A 12 A 2 i = B 12 B 21 + C\ 2 F 2 \ = G\ 2 L 2 \ + Hi 2 H 21 = 1 (3-4) 

as well as 

^12 _ Dl2 Bi2 _ E\2 _ _C\ 2 _ F\2 H\2 _ if 12 _ G*i2 _ _Li 2 

D 2 ± A 2 \ ' £21 -B21 C21 £21 ' K 2 \ H21 G21 L 2 i 

(3.5) 

Given the R-matrix coefficients in Tab. 2 it is easy to convince ourselves that this 
property indeed holds if the undetermined factor R\ 2 obeys the equation R^R^ = 1. 
Curiously, the coefficients in Tab. 2 even satisfy the stronger relation that all three terms 
in (3.5) are actually equal. This amounts to 

A\ 2 D\ 2 = B\ 2 E\ 2 — C12-F12 = B\ 2 K\ 2 — G\ 2 L\ 2 . (3.6) 
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R! 



o g Cl Ui x\ - x 1 
12 q c *U 2 X2 - xf 



q c *U 2 

R o ^ + q -ij ^a _1 727i? Cl t / i ig^xt - (5 - <r X ) - s(xj") 



q 2C 2 +3/2jj2 x ^ _ s ^ x ^ j x ^ _ x j 



- R 12 



1 - (? + q- l )q- 2C *- l U 2 -* X -^\ X l _ X L 

_po / -in ^~ 1 727ig Cl ^i ig -1 4 - (g - g' 1 ) s 04) ~ s ( x f) 

n 12 {q + q ) ^+3/2^2 " x -- s ( x -) ' x^-xt 

a 2 U 2 q C2+1 / 2 (xt - X2) U iq c ^ 2 (xt - x±) 



l-g 2 {q-q~ 1 ) 2 7 2 2 7l 2 

1 ^ 1 



q C 2 +l/2 U2 x - _ x + 

dO Tj. ^ ~ ^2 
-"■12 - + 

72 x 2 - 
q Cl Ui 72 x| -xj" 



R< 



12 



5^2 [/ 2 7l x _ x 



i2? 2 g Cl+1/2 l/i^ 



2 

X o X 1 



:r 2 



Table 2: The coefficients for the fundamental R-matrix of U 9 (fj). 
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Yang— Baxter Equation. Furthermore we have considered the Yang-Baxter equation 

n 12 n 13 n 23 = n 23 n l3 n l2 . (3.7) 

It amounts to around hundred cubic equations among the coefficients A, . . . , L. We have 
confirmed in Mathematica that all relations hold subject to the constraint (2.61). An 
alternative argument |17| for the validity of (3.7) uses the decomposition of the threefold 
tensor product of fundamental multiplets, see (2.35 2.36) 

(0, 0) <g> (0, 0) <g> (0, 0) = {1, 0} © {0, 1}. (3.8) 

Effectively we have to prove the YBE only for one component in each of the result- 
ing multiplets. Representative states for the two multiplets are given by I^VV 1 ) an d 
\%p 1 %p l ip 1 ) . The YBE is trivially satisfied for both and thus it is valid in general. 



Matrix Unitarity. The R-matrix is also a unitary matrix 

(K 12 yn 12 = 1 ® 1. (3.9) 

Using the above unitarity we can rewrite the condition as (Tli 2 y = TZ 2 \ which makes it 
straightforward to read off unitarity conditions for the coefficients 

(A 12 )* = A 2i , (B 12 )* = B 21 , (C 12 )* = F 21 , (G 12 )* = L 2 x, (H12)* = H21, 

(d 12 )* = d 2U (E 12 y = E 21 , (F 12 y = c 21 , (l 12 )* = g 21 , (K 12 y = x 21 . (3.10) 

These conditions are satisfied when the conjugate parameters are given by (2.76) and 
when the phase factor is a pure phase (-R? 2 )* = -^21 = (-^12) 



Crossing Symmetry. Finally, the R-matrix may in principle have crossing symmetry. 
The crossing equation for the fundamental R-matrix takes the standard form 

((T 1 ® i)nf 2 m (c ® i)n 12 = 1 ® 1. (3.11) 

Using the same trick as above for matrix unitarity, we can write the crossing relation in 
terms of the R-matrix coefficients as follows 
a- - A 2 i — B 21 A\ 2 — B\ 2 C\ 2 tt_ ^21 

A 12 — ; Z7~ ' , Z7~ ~ A 21, <J"12 — -^21; : ZT — "21) -"12 



q + q 1 q + q q + q g + g _i 

_ -^21 — E 2 i D\ 2 — E\ 2 F\ 2 C 2 i 

u \2 — : r~ ) : 7— — L) 2\- l -^12 — ^21, — : 7 — ^2i ; -"-12 — — : 7- 

q + q 1 q + q 1 q + q 1 q + q 1 

(3.12) 

All these relations hold simultaneously if the phase factor satisfies the crossing relation 

D D X l X 2 x l b \ x 2 J -l x l x 2 b \ x l J x 2 

II 1 > / 1 . 



H2 ±L 12 



_ ^ g ( X l ) S ( X 2 ) ■s(^i') 3^2 _ g-l s ( x l) S ( X 2) X l S ( X 2 ) 
S ( X l) ~ S ( X t) S { X l) ~ X 2 S { X l) ~ S ( X t) X l ~ S ( X 2~) ' 

which can be written in many (more) alternative ways. This is the quantum-deformed 
analog of the crossing relation obtained in |39|. It would be interesting to find a solution 
of (3.13). 
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Figure 3: Tensor product structure of two fundamental representations along a chain. (1,0) 
contains the state l^ 1 ^ 1 ) while (0, 1) contains the state |V' 1 V' 1 )- The arrows indicate the action 
of the long representation. The long arrows break when 

x<2 — x or — X^ ■ respectively. 

3.3 Special Points 

Let us investigate the behavior of the R- matrix at special points of the parameters xf 2 - 
The analysis is analogous to the undeformed case in |17,40|, here we merely give a brief 
summary. 

Permutation. As usual when the representation parameters coincide, xf — x 2 , the 
R-matrix becomes a graded permutation operator times —1. 

Standard Poles and Zeros. The R-matrix has a pole at x 2 = x\ and a zero at 
x 2 = x i- At these points the shortening condition (2.32) holds, and thus the long 
multiplet in (3.2) splits up as follows (2.31) 

{0,0}-, (1,0)0(0,1). (3.14) 

Note that the action of U 9 (f)) on two fundamental multiplets is defined via the coproduct, 
and typically the resulting long representation is not decomposable, see Fig. 3. In par- 
ticular, at x 2 = x± the representation closes on the short multiplet (1, 0) which contains 
the state 10V 1 ). The residue of the corresponding pole of the R-matrix projects to this 
submultiplet. Conversely at X 2 — X ^ the representation closes on the short multiplet 
(0,1) which contains the state [ip 1 ^ 1 ). The R-matrix projects to this submultiplet at 
this point. This behavior is in fact standard for Yangian and quantum affine algebras, 
but here it appears already at the level of the finite-dimensional algebra U 9 (f)). 

Singlet Pole. Furthermore, there is a pole at x 2 = s(xf) whose residues project to 
the singlet state |1) discussed below. It corresponds to the decomposition 

{0, 0} -> • © adj psu(2|2) © ., (3.15) 

where • represents a singlet. 
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3.4 Singlet State and Quasi- Triangularity 

A singlet state is annihilated by all generators of the Hopf algebra U 3 (f)). Consequently, 
a universal R-matrix in U 9 (f)) <8> U g (f)) would have to act trivially on this state. In a 
quasi-triangular Hopf algebra this also applies to singlet states which are composed from 
non-singlets (fusion). We can use this property to obtain an alternative derivation of the 
crossing relation (3.13) for the phase factor R® 2 as in |13|. 

Singlet State. The singlet component |1) in a two-spin state is defined by A(3)|l) = 
which should hold for all generators Z- Without specifying the spin orientations in |1), 
the constraints A (3) — for the central charges £, ^3, & lead to a relation between the 
representation labels 

xf = s(x 2 ), x± = s(x 2 ). (3.16) 
We make an ansatz which satisfies the constraint for the bosonic generators (Si, £3, t?i, $3 

|i) = - <? +1/2 l</>V> + ^ 1/2 |^V 2 ) - /«7 +1/ >V>- (3.i7) 

The coefficient k is determined by the constraint for the fermionic generators £2 > $2 

a(q Cl Ui - q C2 U 2 ) 

By means of the crossing relation (2.72) for 72 = 71 it simply equals k — 1. 



(3.18) 



Representation of the R-Matrix. The singlet state obeys the curious identity 

^ 13 ^23|ll2X3) = A|l 12 X3) (3.19) 

with a common eigenvalue A for all X 6 {<ft a , tp a }. The eigenvalue takes the form 



3 1 3 2 



A = A 13 (A 23 - E 23 ) = £l3p23 ~ E 23 ) = 4?^^ 

q + q~ l q + q~ l 13 23<? qW 3 g° 3 C/ 3 ^3 - x± x 3 - xf 

(3.20) 

This observation is in agreement with quasi-triangularity of the Hopf algebra U 9 (f)): If 
we view the two-particle state as transforming under one representation (the singlet) 
then the product 7£i 3 7?.23 must equal the representation 7?.i3 of the R-matrix on a singlet 
and a fundamental (fusion). However, singlet representations of the R-matrix are trivial 
and thus the eigenvalue must be A = 1. This condition leads to a constraint on the 
overall phase factor 

rf rf rf ' q ( rf ~\ 

R 12 R I2 - q — _ . (3.21) 



It is the same as the crossing relation obtained earlier in (3.13). Note that we have 
mapped the representations according to (1, 2, 3) — > (1, 1, 2). The map (1, 2, 3) — > (1, 1, 2) 
leads to a contradictory result. A consistent result is achieved by (1, 2, 3) — > (1, 1, 2) for 
which we obtain 

pO pO _ _ S ( X 1 ) ~ X 2 s ( x l ) ~ s ( x 2 ) (a nn\ 
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Note that B2 ^ R\ 2) more precisely there is a double crossing relation which is the 
quantum-deformed analog of the relation solved in 41 

-^12 _ x i ~ x t s ( x i ) ~ x 2 x t ~ s ( x t) s ( x t) ~ ■ s ( a; 2 ) ^ 



4 Diagonalizing the R-matrix 

The R-matrix can be used to determine the eigenstates for certain spin chains. In 
particular, it is required to write down quantization conditions for eigenstates of finite 
closed or open spin chains. However, as the R-matrix is a matrix, the resulting equations 
would be matrix equations and rather hard to handle. Instead of dealing with matrix 
equations we can first "diagonalize" the R-matrix by introducing a suitable vacuum state 
and excitations. This is done by means of the nested coordinate Bethe ansatz |42|, and 
the procedure follows along the lines of |17| for the undeformed setup. Alternatively, one 
may perform a nested algebraic Bethe ansatz similar to the one in 1 5 , 43 1 which should 
lead to the same set of Bethe equations. 



4.1 Vacuum and Propagation 

We would like to find suitable eigenstates of a chain of K fundamental representations 
with labels x^, k = 1,...,K. We define the level-II vacuum state |0) to consist only 
of 4> lJ s. The top-level R-matrix TZ\. is a product of pairwise R- matrices = T^k,j 
representing a permutation 7r G Sk of the K sites. It multiplies the vacuum state by a 
phase factor 

KW l = R\W l - (4.i) 

The total phase factor is a product over pairwise phase factors 

K= IT Rl ' l ^k,Xj), R l >\x k , Xj ) = A(x k , Xj ) = R°(x k , ^)^7r X - ~ X l ■ ( 4 - 2 ) 

(fcj)S7T J J 

The level-II vacuum has two flavors of elementary excitations. An excitation replaces 
one of the <ft l spins by a tf) a spin 

lf)! = EM»)l^-f-ff- ( 4 - 3 ) 

k=l 

The wave function \Pn,k{y) is parametrized by a rapidity y and it depends locally on the 
parameters xf of the spin representations 

fc-i 

= f(y, x*{k)) II i?n ' 1 ^' M?))- ( 4 - 4 ) 

i=i 

The permutation tt defines the ordering of the top-level spins along the chain. 
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In the nested Bethe ansatz the wave function must be compatible with the action of 
the R-matrix 

K\rf = Rl\r) l l (4.5) 

It suffices to consider a chain with two sites and the permutation 7r interchanges them. 
The solution to the compatibility condition reads 

R ll \y,x k ) = q c ^u k y -—^ , f{y , Xk) = . (4.6) 

Clearly the rapidity y could be replaced by any function of y; we chose it such that the 
functions R ll,l {y,Xk) and f(y,x k ) are rational. Furthermore f(y,Xk) could be multiplied 
by any function of y. 

4.2 Scattering 

We now have to solve the analog of the compatibility condition (4.5) for a two-excitation 
state. We make the following ansatz 



vr(fc) 7r(j) 

|VV)?= E *,jA!h)*:,,(!)i)^ >:->... <y) x (4.7) 

k<j=l 

which solves the compatibility condition except when the two excitations are neighbors. 
This state can mix with a state where one spin 1 is replaced by (f> 2 

\<PX = E #«*(yi)#*Mf(vuy2, ^(k)M ■ ■ ■ <P 2 ■ ■ ■ <f> 1 ) 1 - (4.8) 

k=l 

The spin <p 2 should be interpreted as the overlap of two fermionic excitations. Thus a 
generic two-excitation eigenstate reads 

i#>; = + c"*v>" + K^ir^l (4.9) 

The matrix C a ^ has the non-zero elements C 12 = q^ 1 ^ 2 and C 21 = —q +1 ^ 2 . The operator 
V\\ is a (graded) permutation which interchanges the excitations along with their rapidi- 
ties 2/1, 2/2- Finally the level-II R-matrix TZ 11 must be invariant under U 9 (su(2)) which 
restricts its form to 

K ll \^ 1 ) 11 = M l2 \^ l ) l \ 

-olll il / 2\II — ^Vl2 | ,i , 2\ II , + q~ l N\2 1,2/ 1\II 

q + q 1 q + q 1 

-0111/2/1x11 <T l M 12 + qNi2, ,i 2 \ii , M i2 - ^12, ,2 

U |V V ) = q~^q=i + + i W 1> ) > 

^ n |VV) n = M 12 iVV) n (4-10) 
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Imposing the analog of the compatibility condition (4.5) yields a set of equations 
from which the unknown functions M^yi, 1/2), -^12(2/1, 2/2), f(yi, IJ2, x k ) can be extracted. 
The level-II R-matrix 1Z ' can easily be determined by considering any of the triplet of 
states |V>V) n , q 1/2 \^V) U + <T 1/2 |^V> n or \ip 2 ^ 2 ) 11 . It follows that 

M( yi ,y 2 ) = l. (4.11) 

To determine the other functions we pick the singlet state g~ 1//2 |'?/> 1 -?/> 2 ) 11 — q 1 ^ 2 ^ 2 ^ 1 ) 11 . 
A lengthy calculation shows that 

v qu(y 2 ) - q~ l u(yi) + ig~ l 
N m, V2) = =1 7 s ^—r r-zi , 

fr« 7, _ + g~V /2 xt{xt-xp yi - a/2 

2 ^7fc ig 1 + {q-q x )x k q l u(y 2 ) - qu{yi) - ig 

„r„ \ a /„ \ VlV2 ig ~ l + ^-^ r ) x k i^ 1 + (q- q~ l )xj 
< x k % g vi -\q-q ) w yz - (? - 9 ) 

is the solution to the compatibility condition if u(y) is the function, see (2.73), 

u(y) = y + s(y), s [y = — 3—. 4.13 

w y-{q-q ) 

Note that the level-II R-matrix indeed solves the YBE and u(y) is the rapidity variable 
for this Ug (bu(2)) R-matrix. The function u(y) is invariant under taking the antipode 

u(y) = u(s(y)), (4.14) 

because the map s(y) is an involution. 

4.3 Final Level 

The final-level R-matrix has the standard form for U 9 (su(2)) symmetry. The diagonal- 
ization by means of the nested Bethe ansatz leads to the following phase factors 



q C] U 3 H ~ x j 
q Ck U k x7 - x+ ' 



R ll '\y v x k )=q c ^/ 2 U k ^^, 

R u ' ll ( yj ,y k ) = M 12 ( yj ,y k ) = 1, 
{ w j,yk) — q 



q l u(y k ) - Wj - \g~ A 



1 



D in in/ \ q w k - qWj - ^{q + q x )g 1 

R m ' m { Wj ,w k ) = 2 . (4.15) 

qw k -q 1 w j + ^(q + q l )g 1 

Here w k are the level-III rapidities. For completeness we have included the previously 
derived phase factors from lower levels. 
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4.4 Bet he Equations 



Let us now consider a closed spin chain with K sites. The wave function for excitations 
must be periodic in order to define suitable eigenstates. Periodicity is imposed by means 
of the Bethe equations which use the elements of the diagonalized R-matrix. Generically 
for a nested Bethe ansatz with two levels they can be written as follows 



K N M 



j=i j=i j=i 

K N M 

1 = ]\R 1 ' n \x v w k )]lR u ' n \y 3 ,w k )]lR m ' n \w J ,y k ), (4.16) 

j=i j=i j=i 

Thus for the system in question they take the form 

K + M ( n, i _i 



l = ff q -c i -i/2 U -x yk- x o tt q -i qujyk) - Wj + 
jJi 3 Vk - x] fJi q~ l u(yk) - Wj - fir 1 ' 

1 = tt Wk ~ q~ lu (yj) + 1^ 1 tt q' lnj k - q w j - Kg + g' 1 )^ 1 



(4.17) 



4.5 Dualization 

At the very beginning of the nested Bethe ansatz one can in fact choose between two 
different vacua, composed from only bosons 1 's as above (alternatively 2 's) or composed 
from only fermions ip l, s (alternatively -^ 2 's). In both cases the NBA will proceed in a 
very similar fashion, but lead to two different but equivalent sets of Bethe equations. 

Instead of performing the alternative NBA, we shall derive the alternative Bethe 
equations by means of dualization. The dualization of our Bethe equations for the 
undeformed case was performed in |18| (which is equivalent to the dualization of the 
Lieb-Wu equations |3| in |44|). Here the procedure is qualitatively the same, but requires 
special care due to various insertions of q. Let us outline it in the following: 

The first Bethe equation in (4.17) can be viewed as an algebraic equation P(y k ) = 
with the polynomial 

K M 

P{y) = II q Cj+1/2 u, (y - xj) J] q {ig-'y -( g - q- 1 )) (<rV y) - Wj - \ g - 1 ) 

i=i j=i 

K M 

- n (v - x t) II i^'y -(9- q' 1 )) {qu{y) - ^ + t^r 1 ) • (4.is) 
3=1 3=1 

Note that the common factor (ig~ x y — (q — q^ 1 )) was introduced in order to cancel the 
denominator of s(y) in u(y). This polynomial has degree K + 2M, i.e. of its roots 
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are the yu and the remaining N = K + 2M — N roots will be denoted by yj.- We can 
construct a constant function F(y) by dividing by all root monomials 



N . N / K 



nv) = p(y) n zhf. n ^j. = iv>- i ) M n ^ +i/2 ^ - • ^ 

i=1 f % \i=l / 

Next let us define as a function of implicitly through the relation 

w k = q' 1 u(x+) - \g~ l = qu(x~) + \g~ x . (4.20) 

The x^j obey the x ± -constraint (2.61) which can be written in a particularly convenient 
form using the map u(y) (4.13) 

q~ l u(x + ) - qu(x~) = - . (4.21) 

9 

Now any ratio of the constant function F(y) evaluated for two different values of the 
parameter equals 1. There are four useful points y = x^,s(x^) where one of the two 
terms in the polynomial P(y) drops out by means of (4.20). At these points we obtain 
the following identity which holds by virtue of the first Bethe equation in (4.17) 

1 = F(x+)F{s(x+)) = Y^ q - 2C ,-i U7 2 x t ~ x t s ^ - x t 
F(x w )F( k s( k x w )) ^ 3 x w — Xj s(x w ) — Xj 

TT x w ~ Vj s ( x w) - Vj yr x- - yj s(x~) - yj 
jJl x t - Vj s ( x w) - Vj jJl x t, - Vj s ( x t) - Vj 

- 2 w~ lx w - (? - 'T 1 ) ig~ x s(x+) -(q- q' 1 ) 



q 



w lx w - (q - q x ) w ls ( x w) -(q-q x ) 



^ Vg- 1 ^^-) - Wj - §#- ] / 



For the first line we use the identity 



J^^_ ^ = g 2C J+2f/ 2 (4 23) 

^ - S i S K ) - X j 

and obtain simply g K . The second line is simplified by means of the identity 

x + — y s(x + ) — y u(x + ) — u(y) 
x~ — y s(x~) — y u(x~) — u(y) 

and for the third line the identity 



(4.24) 



ig './'• - (q-q l j o r , r2 



ig~ 1 x — (q — q x ) 
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q 2C U 2 (4.25) 



leads to q 2M . In the last line we substitute the definition (4.20) of x^. Altogether this 
yields the identity 

1 = TT -i w k - qu(yj) - \g~ x yj _ 1 w k - qujyj) - \g~ x 

7=i 9 w k- Q'Mvj) + y- 1 jLi 9 Wk - <?~m%') + y^ 1 

jrf qw k - q- l Wj + |(g + q- l )g- 1 \ 2 ^ 



Two of these terms coincide with terms in the second Bethe in (4.17). We multiply the 
Bethe equation and the identity and obtain the dual Bethe equations for y k and w k 

i = fr q Ci+v*u- ^ k - x ~ j fr q q ~ lu ^ k) - Wj - ^ 

1 = TT g -i w * ~ g^fe) ~ l9~ l yr qw k - q~ l Wj + |(g + g" 1 )^ 1 ^ ^ 

\}[ w k - q^uiyj) + ^g- 1 jJ^q^Wk- qwj - ^(q + q-^g- 1 ' 

These Bethe equations do not contain the scattering term R}' l (xj,x k ) which should 
be dualized as well for completeness. In fact it can easily be obtained from the scattering 
in the dual level-II vacuum composed from ^ s: it equals —D(xj,x k ). An alternative 
procedure is to consider the x^ as dynamical degrees of freedom as in 1 18 1. In other 
words they obey some Bethe equation of the form 



K 


N 


M 






>n 


3=1 


3=1 


3=1 


3^k 






K 

JJi2°(xj,x fc ) 

3=1 


q c Wj xt -xj 


N 

n 

3=1 


q Ck U k x k - xj 


3+k 







q C k +l/2 Uk X k_Jh (4>2g) 

x l ~ Vj 



We can now derive another identity from the constant function F(y) 

n/ -\ N i N + _ K -4 

x = F(x k ) = _ TT x k ~ V_i TT x k ~ Vj TT q -C--l/2 U7 l X k ~ x j 

F ( x i) f = i x k - Vj x k - Vj 3 x i~ x j 



q 



r- 

M M I —\ i i -1 



ig 1 x+-(q-q x )\ qu(x k ) - Wj + \g 



Lln-U,{rr+\-,n.-Ln-l- K > 



W~ l x k - (q - q' 1 ) J fJi q- l u{xl) - Wj - \g~ x 
Using the identities (4.21.4.25) we find 



N + N + ~ K C k TT - — + 

i = n^-^^n^- 1 "^ 1 ^! n y£ S , ( 4. 3 o) 
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which yields when multiplied to (4.28) the dual Bethe equation 

k n + ~ 

t. ^ TT n -c k -i/2 TT -i x k - Vj 



n ^ n g-^- 1 / 2 ^- 1 ^^ . (4.3i) 
j= i i= i % 



In conclusion the dual phase factors read 

(xj^Xj^j D^Xj^Xj^j R (xj,x^) ; 

x k ~Vj 

K (w i; y fe J - q — — , 

5111,111/ \ gwfc-g~ 1 ^ J + Kg + g" 1 )^ 1 

i? K-,w fe ) = — — - . (4.32) 

q L w k - qwj - ^{q + q L )g 1 

5 Quantum Deformation of the Hubbard Model 

In this section we derive an integrable homogeneous nearest-neighbor Hamiltonian from 
the U g (f)) R-matrix in Tab. 1. Its Hilbert space is the same as for the one-dimensional 
Hubbard model and we will show that our Hamiltonian is a three-parameter deforma- 
tion of the Hubbard Hamiltonian. We then discuss the relation to the integrable two- 
parameter deformations of the Hubbard model proposed by Alcaraz and Bariev |34| and 
introduce additional sets of two-parameter deformations. 

5.1 Hamiltonian 

In the following we derive a nearest-neighbor spin chain Hamiltonian from the U g (f)) 
R-matrix in Tab. 1. 

Integrable Hamiltonian. For an integrable spin chain based on an R-matrix there is 
a standard procedure to obtain a homogeneous Hamiltonian 

L 

k=l 

The pairwise interaction H.12 is the following logarithmic derivative of the R-matrix 



. (x+ - s(x+)) (x~ - s(x~)) fdu*Y 1/2 ^ , d „ 
q 1 x + s(x + ) \du J dui 



(5.2) 



The spectral parameters Uk are defined as functions of the parameters x k , see (4.13), 

u k = q~ 1 u{x+)-^ = qu{x k ) + ^. (5.3) 
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n 



121 



Hl 2 |0V> 



^12 

H-X2 
7~tl2 



n 



12 



Uvi\rv) 



^+V« + 

g + g _1 



c 



g + g 1 

Ala^ 2 



g + g 1 



i^ + 9 1 ^ 4 + ^ l ^ 2 ^i\ 



g + g 



g + g x ' g + g 

g + g 



) 

£>|V>V> 
qD + q^E 



IV>V 2 > + V> + " 



-1^ 



g + g- 



D-E. ^ , ?^ + 9 B^ V) 



g + g 



g + g 



-i 



g + g 1 

g + g _1 



G^V) + > 



g + g" 



;li2 



2 ) + 



g + g 1 
qF 

q + g _1 



i2il\ 



12 /h 



Table 3: General form of the nearest-neighbor U 9 (su(2) x su(2)) spin chain Hamiltonian 



Note that u k is not real 



n't 



(g + g 1 )u k + |g x (g-g x ) du£ 
(g + g -1 ) - 2ig(q - q~ x )u k ' ciw fc 



(g 1 x+s(x+)) 2 =(gx fc s(x fe )) 2 , (5.4) 



and therefore we have included a compensating prefactor (du* / du)^ 1 ^ 2 = q^ 1 x^s(x^) to 
make the Hamiltonian hermitian. The additional real prefactor is meant to simplify the 
resulting expressions. 

The resulting Hamiltonian clearly has U g (su(2) x su(2)) symmetry 6 and the general 
form with such symmetry is listed in Tab. 3, Just like the R-matrix in Tab. 1 is has ten 
independent coefficients X = A,B, C, D, E, F, G, H, K, L. We have arranged them such 
that they are directly related to the corresponding coefficients Xu of the R-matrix in 
Tab. 2 according to 



X = i{x + — s(x + )) (x — s(x )) 



dX 12 



du\ 



, _,_ _,_ , , _ / dx 0X\> dx dXi2 



(5.5) 



6 It does not necessarily have U g (P)) symmetry due to the u-dependence of the coefficients a, b, c, d in 
the representation of the fermionic generators. In fact there is no U 9 (t))-invariant operator apart from 
the identity for generic q. 
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_1_ (q C U + q-CU-^jqCU- 1 + q- C U) 
Ag (q c U - q-CU-^iqCJJ- 1 - q~ c U) 



g (q c U - q- c U- l )(q c U- 1 - q-°U) 
C = F=(g + g- 1 ) v /l-(g-g-i)V 

q C+l/2jj-l _ q -C-l/2jJ-l _ q C-l/2 U + q -C+l/2jJ 

G= 9(q- q~ 1 )(q c U - q- c U~ l )(q c U- 1 - q~ c U) 

q C+l/2jj _ q -C-l/2jj _ q C-l/2jj-l + q ~C+l/2jj-l 

L = g(q - T^fU - q-CU^iqCU- 1 - q~ c U) 
H = K = 



Table 4: The coefficients for the nearest-neighbor Hamiltonian. 



Note the following two useful identities in evaluating the derivatives 



dx ± x + - x~ o ±c ' ±1 t/ ±1 



du x ± — s(x ± ) q c U — q C U 



-l ' 



4- log(Uq c ) = - ( - q , q ^T—^\ ■ (5.6) 

du 6V H ' 2 \x+-s(x+) x--s{x-)J y ' 

The coefficients X still depend on the parameters a, 7^ and the phase factor R® 2 . 
The latter two can furthermore depend non-trivially on x^ which would influence the 
Hamiltonian explicitly or implicitly through the derivatives. For definiteness we set the 
phase factor to 

a different phase factor would induce an overall shift of the spectrum which we shall 
incorporate explicitly later. The most suitable expression for 7^ is given in (2.65) and 
we shall set the global parameter a to unity, a — 1. 

Now we are ready to obtain an explicit expression for the nearest-neighbor Hamil- 
tonian. Taking into account that for xf = xf the R-matrix becomes minus the graded 
permutation operator one can show that the Hamiltonian has the coefficients listed in 
Tab.S 



Integrability Constraints. The coefficients obey certain relations: There are two 
linear relations and two quadratic relations 

A + D = B + E = H + K, 
[A-Bf + CF ={q + q- 1 ) 2 GL, 
2{A-B){A-D) = {q + q- 1 )(G 2 + {q + q- 1 )GL + L 2 ). (5.8) 
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Note that the two linear relations can be derived from the identity (3.6); the origin of 
the two quadratic relations remains unclear. These equations are invariant under four 
trivial transformations: (i) A rescaling of all coefficients by a common factor, (ii) A 
shift of the two-site Hamiltonian by the two-site identity operator A(l), taking the form 
in Tab. 3 with 

A = B = D = E = H = K = 1, C = F = G = L = 0. (5.9) 

(m) A reciprocal rescaling of C and F corresponding to a different rescaling of fermions 
\ip a ) w.r.t. bosons \<fi a ). (iv) An opposite shift of H and K which has no impact on the 
spectrum. Altogether the ten coefficients together with the parameter q subject to the 
four constraints (5.8) and four identifications (i-iv) define a three-parameter family of 
models. This is the same number of degrees of freedom as for the fundamental R-matrix 
given in terms of q,g,x ± subject to the one constraint (2.61). Therefore the constraints 
(5.8) are expected to be sufficient to ensure integrability of the Hamiltonian in Tab. 3, 



Hermiticity. There are ten independent complex coefficients in the Hamiltonian in 
Tab. 3, The hermiticity condition 7i\ 2 = imposes certain reality conditions on the 
coefficients 

C* = F } G*=L, A,B,D,E,H,K <E R. (5.10) 
Furthermore we have to require that 7 

qER. (5.11) 
Due to our choice (2.65) for 7*., we have that C = F and thus C,F must both be real. 

Generalizations. Before performing an explicit comparison of the models in question 
we would like to introduce certain integrable generalizations of the Hamiltonian which 
change the spectrum in a controllable fashion. Thus given a Hamiltonian 7i 12 as in Tab. 3 
one can transform it as follows 

H'u = a THuT- 1 + iaiACQi) + a 2 A(l) + ±a 3 A(S) 3 ) (5.12) 
+ §&i(#i ® 1 - 1 ® £1) + 62(^1-61 ® 1 - 1 ® Mi) + lh{Sj 3 g> 1 - 1 <g> Sj 3 ). 

Here a&, bk are arbitrary constants. Indeed a$ is an overall multiplier of the operator, a 2 
is an overall shift and a±, correspond to a shift of the energy eigenvalues proportional 
to the $Ji,$)3 Cartan generator eigenvalues. The terms multiplied by bk vanish after 
summation over the whole spin chain with periodic boundary conditions (5.1). 

The similarity transformation induced by T is the following Reshetikhin twist |45| of 
the integrable structure 

A' (3) = TA(3)T-\ K' = V{T)nT- 1 . (5.13) 

7 Usually in quantum-deformed spin chains the case of q on the unit circle also leads to hermitian 
Hamiltonians. Here we have not been managed to establish hermiticity in this case. 
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It is not hard to verify that the twist preserves the cocommutativity property. We shall 
consider a twist T constructed from the identity and Cartan generators 

T = exp(i/i(l <g> $i) + §/ 2 (£i <g> io 3 - £3 ® fli) + i/ 3 (l ® «s)) , (5.14) 
which can be generalized consistently to arbitrarily many sites 

(K k k > 

i=i j<fc=i i=i 

(5.15) 

The coefficients /& are arbitrary parameters. Hermiticity requires them to be real. 
5.2 Bet he Equations and Spectrum 

The spectrum of the above Hamiltonian on a closed homogeneous spin chain is deter- 
mined by Bethe equations. Here we specify the Bethe equations and energy relations. 

Energies. We use the Bethe ansatz based on a ferromagnetic vacuum consisting of K 
spins 1 as in (4.1). We assume there are N main excitations (magnons) with momenta 
Pk which turn a <p l into a ip 1 . Finally, there are M auxiliary excitations with rapidities 
Wk which turn a ip l into a ip 2 . We have shown in Sec. 4 that the Bethe equations for 
this system are given by (4.17) where we set all representation parameters to be equal 
x fc = x± f° r a homogeneous chain 

-G-\l2jj-\ Vk -x + \ -pj- _ x qu{y k ) - wj + ~ 2 g 



1 _ TT w fc -g l u(y.j) + lg 1 ji q 1 w k - qwj - |(g + q l )g 1 

J = i Wk~ qu(yj) - \g- 1 qw k - q^Wj + |(g + g" 1 )^ 1 



(5.16) 



First of all, a relation between the magnon momenta and the magnon rapidities 
y k has to be established. In the Bethe equations (4.16) the term B} ,ll (x,yk) serves the 
purpose of e~ ipk , i.e. the rapidity relation is 

N 

e i Pk = R^\ yk)X ) = q c+V*u y±^^ e iP = TT e i Pk . ( 5.i 7 ) 

y k — x + A A 

yfc j=1 

Here the total spin chain momentum is given by P. The energy for a solution of the 
Bethe equations is given as the sum of vacuum energy and magnon energies (here K 
refers to the length of the chain) 

N 

E = E K + J2E(yj). (5.18) 
i=l 
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The vacuum energy density Eq and magnon dispersion relation E(y k ) follow readily from 
the expression for the Hamiltonian in Tab. 3 (here K refers to the coefficient K of the 
Hamiltonian) 

E = A, E(y k ) = H + K-2A + Ge ipk +Le~ ip \ (5.19) 



where e ipk is defined in (5.17) as a function of y k . The energy relation can also be 
obtained formally from the diagonalized elements of the R-matrix in analogy to (5.17) 



Eq = i(x + — s(x + ))(x — s(x )) — loglth l (xo, x) 

du 



E(y k ) = i(x + - s(x + ))(x- - s(x~)) — \ogB} l >\y k) x). (5.20) 

du 

This leads to the same as the above expressions. 

Transformation. The above generalization of the Hamiltonian (5.12) requires certain 
modifications of the just derived expressions. Most importantly, the Bethe equations 
receive further phase factors due to the twist induced by the f k 8 

iih-h-h) q -c-in u -i yk-x + \ K ^r c2ih ^ l qujyJ-Wj + ig- 1 

Vk-x-J q- x u{y k ) - Wj - \g~ x ' 

fJl w k - qu(yj) - \g- x jJ^qWk-q^Wj + liq + q-^g- 1 ' 

3+k 

These phase factors therefore enter the momentum relation e tpk = R Il,1 (y k ,x) as follows 

e i Pk = e i(/i+/2-/ 3 ) g c + i/ 2[/ VkZ^L . (5.22) 

Vk ~ x+ 

Due to the logarithmic derivatives in (5.20) the expression for the energy is unaffected 
by the phase factors. It however receives contributions from the rescaling and shifts in 
(5.12) ' 

N 

E' = (aoEo - ai + a 2 )K + 2a 3 M + ( a oE(yj) + oi - o 3 ) . (5.23) 

i=i 

Dual Picture. The procedure to obtain the spectrum in the dual picture discussed in 
Sec. 4.5 is the same. Here the ferromagnetic vacuum is based on ^'s, the magnons flip 
them to 1 's and auxiliary excitations flip the latter to 2 's. The twisted version of the 



^Clearly the twist can also be applied to the inhomogcneous equations where x k ^ x, ■ 
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dual Bethe equations in (4.27) have the following insertions of phase factors 

-\ K M -l (~ \ _ _ i _i 

Jih+h-Mf+ifiU yj^l^) TT e -*A q q y>) w * & (5.24) 

jLi Wfc ~~ + \9~ l fJl q~ l w k - qwj - |(g + T 1 )^ 1 ' 

To obtain the correct momentum relation we have to remember that the vacuum as well 
as the magnons are fermionic excitations. Therefore there is an implicit factor of —1 in 
the momentum relation e l ^ k = —R ' (yk,x) 



N 



e ip k _ _ e i{f3-h-h) n -C-\/2jj-l Vk x 



q-c-i/tjj-i yj^_^_ e w = (-l) K+N - 1 T\e i ^. (5.25) 
y k -x a± 

Note that the definition of the total momentum P via the cyclic shift operator requires 
taking into account the statistics of the vacuum sites and excitations. The energy relation 
uses the above dispersion relation E(y) up to a sign 

TV 

E' = (-a E -a 3 + a 2 )K + 2a x M + ^ {- a o E (Vj) - ai + a 3 ) . (5.26) 

3=1 

The equality of the total momentum P and total energy E' in both pictures makes use 
of the identity (4.30) which guarantees for all eigenstates that 

N N N N 

j=i j=i j=i j=i 

5.3 Electronic Oscillator Notation 

We have established the Hamiltonian in terms of a matrix acting on a spin chain based 
on four-dimensional vector spaces. Let us rewrite the four-dimensional vector space in 
terms of fermionic "electron" creation and annihilation operators and c a as it is usually 
done in the Hubbard model literature. 

The map between vector states and electronic states reads 9 

= |o>, \<f> 2 ) = Kc\ci\o), IV 1 ) = 41°), |^ 2 ) = 4|o). (5.28) 

where creation and annihilation operators satisfy the following algebra 

K, 4} = Sap, {c a , c p } = {cl 4} = 0. (5.29) 



9 More generally, we may introduce different normalization factors for the four states. By rescaling 
the state |o) and the generators c\,c\ these can all be absorbed into a single variable n. 
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We define the number operator by 

n a = c ] a c a . (5.30) 

Note that n 2 a = n a , hence n a is a projector onto the subspace spanned by \ip a ), \4> 2 )- The 
following projectors can be specified 

\<f> 1 )(<f> 1 \ = {l-n 1 ){l-n 2 ), 
|0 2 )(0 2 | = n x n 2 , 
\^ 1 )^ l \=n l {l-n 2 ) ) 

I^X^I =n 2 (l-m). (5.31) 

From these projectors many useful combinations can be assembled straightforwardly. 
The Cartan generators can be rewritten in terms of number operators as follows 

S)i = ni + n 2 - 1, S)3 = n 2 -ni. (5.32) 

For dealing with spin chains we introduce multi-site oscillators by adding a site index 
j, k. The oscillator algebra becomes 

{ c a,P c l,k} = ^fc<W, {c aJ , c Pjk } = {cl j7 4 J = 0. (5.33) 

By means of the projectors specified above we can spell out the Hamiltonian of our model 
in oscillator notation. Thus for the operator from Tab. 3 we have the expression listed 
in Tab. 5, 

5.4 Connection with the One-Dimensional Hubbard Model 

In two special limits the spin chain with undeformed f) symmetry becomes equivalent to 
the one-dimensional Hubbard model |17|. This conclusion was drawn first by comparison 
of the Bethe equations of both models and second by comparison of their R-matrices. 
Here we show explicitly how the Hubbard model can be embedded into the f) spin chain 
on the level of Hamiltonians as a preparation for the more complicated comparison with 
the Alcaraz-Bariev models. 

The problem in the comparison of Hamiltonians consists in the fact that the original 
Hubbard model Hamiltonian exhibits manifest su(2) symmetry while our Hamiltonian 
has manifest U g (su(2) x su(2)) symmetry. We therefore need to recover an additional 
manifest su(2) symmetry and set q — 1 before the Hamiltonians can be compared. 

Hubbard Hamiltonian. The standard Hubbard model Hamiltonian based on the 
electronic oscillator algebra (5.33) reads 

n ff = E ( c L' c ^ + i*^) + Un ^2,r (5.34) 
0=1,2 
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1 ( C l,j C 2,j C 2,k C l,k + C l,k C 2,k C 2,j C l,j) 



q + q 
D-E 



f K C l,j C 2,k C 2,j C l,k + C l,k C 2,j C 2,k C l,j 



T J 



q + q- 

+ q + q -i C L- C i,fc(9~ 1/C ~ lc '( 1 ^ n 2,j)n 2 ,k - qnFn 2 j{l - n 2 , fc )) 

+ q + q -i c h c iA K ~ lc ( l ~ n hj) n i,k - KFmj(i - m ifc )) 

+ g +^-1 c lfc c ij(9~ 1/cF ( 1 ~ n 2j)n 2 , fc - qK^Cnvjil - n 2lfc )) 
+ ~ 4,fc c 2,j (^(1 - nij)wi )fc - «; _1 Cnij(l - m >fc )) 

+ ctjCi,*^ 1 - n 2,i)(l - n 2 ,k) - Ln 2 jn 2 ,k) 
+ c \,j c 2,k{G{\ - n hj )(l - n ltk ) - Lnijrai,*) 
+ c U c i,j - "2j)(l - "2,*) - Gn 2 jn 2 , k ) 

+ A + (K - A)(nij + n 2 j) + (H- A)(n ljk + n 2>k ) 
+ (A + D - H - K)(n 1:j n ltk + n 2jj n 2jk ) 



+ 
+ 
+ 
+ 
+ 
+ 



A-2H+ qA + q B } n hk n 2 , k + ( A - 2K + - A + qB ^j n 1:j n 2J 



A-H -K + 



q + q 1 

qD + q- x E 



q + q~ 



q + q~ l 

, , T ^ q^D + qE" 

ni :j n 2:k + [A- H - K + q + q _ 1 J n 2J n hk 



-A- D + 2H + 2K 
-A- D + 2H + 2K 
-A- D + 2H + 2K 
-A- D + 2H + 2K 



qA + q^B 


q~ 1 D + 


qE 


q + q" 1 


q + q- 


-l 


qA + q- x B 


qD + q~ 


- X E 


q + q^ 1 


q + q- 


-l 


q~ l A + qB 


qD + q 


~ l E 


q + q- 1 


q + q- 


-l 


q~ l A + qB 


q~ l D + qE 


q + q- 1 


q + q- 


-l 



n 2 ,jn hk n 2 , k 
nijn hk n 2yk 
ni,jn 2J n 2tk 
J n hj n 2J n ltk 



+ (3A + B + 3D + E - AH - AK) n ltj n 2J n hk n 2tk , 
Table 5: Oscillator form of the U g (su(2) x su(2)) spin chain Hamiltonian. 
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It has one manifest su(2) symmetry and a twisted su(2) symmetry. In order to match with 
our Hamiltonian in Tab. [3|4| we should first make the second su(2) symmetry manifest 
by applying a transformation (5.12) 

Wj,k = T j,k H lt" 1 ' 3 ~k + \ u ( l ~ n U ~ n i,k - n 2,j - n 2 , k ~ 2nij-n 2 j + 2n lifc n 2iA; ) (5.35) 
with the twist 

T jtk = exp(|7r(n lifc + n 2 , fe - 1)) = exp(§7r#i )& ). (5.36) 
The resulting Hamiltonian reads 

w i* = S (~ ic L' c «,fc + ic a,fc c «j) + 1^X1 ( n V n V + (1 - ni^)(l - n 2/ ) - |). (5.37) 

a=l,2 t=k,j 

It has manifest U g (5u(2) x su(2)) invariance with g = 1 and takes the form in Tab. 3 
with 

A' = B' = —D' = —E' = \U, G' = —L' = i, C, F' = -2(zk) ±1 , H' = K' = 0. 

(5.38) 

Comparison of Hamiltonians. It was shown in |17| that the Hubbard model corre- 
sponds to the values (0, oo) or (oo,0) of the parameters (x + ,x~). Moreover we have to 
set q = 1 for the undeformed setup. 

In fact the Hamiltonian is singular at these parameter values and there are several 
ways in which the singular limit can be taken. For instance, we can first set q = 1 and 
then take the limit 

(x + ,x-) -f (0,oo). (5.39) 
Alternatively one can take the limits 

{x + ,x~) -> (0,s(0)) or (x + ,x") -> (s(oo),oo) (5.40) 

and only afterwards q — > 1. We rescale the Hamiltonian to match the above G' 

K, [, Hp,.. (5.41) 

This leads to, see (5.9) 

A' = B' = —D' = —E' = — , G' = -U = i, C = F' = -2, H' = K' = 0. (5.42) 

which agrees exactly with the twisted Hubbard Hamiltonian (5.38) upon identification 
of k = — i and 

U = - . (5.43) 
9 

5.5 Connection with the Alcaraz— Bariev Models 

Here we study the relation between our model and an integrable spin chain proposed by 
Alcaraz and Bariev |34| (AB model). 



37 



The Alcaraz— Bariev Hamiltonian. Alcaraz and Bariev |34| proposed a spin chain 
model based on the electronic states introduced in Sec. 15.31 The Hamiltonian is the 
following general deformation of the Hubbard Hamiltonian (5.34) 

ftf* = ( C lJ C l,fc + C \,k C l,j)( l + ^l n 2,j + t 12 U 2)k + t'^j^k) 

+ ( c h c 2,k + 4,it c 2,j)(l + foinij + t 22 n 1)k + t 2 n ltj n 1)k ) 

+ J( C l,j C 2,k C 2,j C l,k + C l,k C 2,j C 2,k C l,j) + ^p( c l ,j c 2,j C 2,k C l,k + C l,j C 2,j C 2,k C l,k) 

+ Vunijn^k + V 12 n ltj n 2 ,k + V 21 n 2)j n 1M + V 22 n 2}j n 2jk + Un ld n 2d 

I rr(l) . T/ (2) T/ (3) . rr(4) 

+ V3 1 'n 2tj n 1>k n 2>k + V3 'n hj n^ k n 2 , k + 'n ld n 2tj n 2tk + 'n ltj n 2tj n ltk 

+ V 4 n ld n 2 jn lik n 2:k , (5.44) 

where 

tn = £4 — 1, £12 = £3 — 1, t x = £5 — £3 — £4 + 1, , _•, 

/ (5.45) 
hi — t\ ~ 1) t 2 2 = ^2 — 1, t 2 = £5 — ti — ^2 + 1- 

It was found to be integrable in four cases which we will denote by A and B . The 
coefficients are related as follows in the case A ± 

ti = et 2 = t 3 = e^4 = sini?, i 5 = e = ±1, 
J = -et p = -\eU = V 12 e 2n = V 21 e- 2r > = costf, 

= ^22 = ^ 3 (1) = 1/ 3 (2) = 1/ 3 (3) = 1/ 3 (4) = 14 = 0, (5.46) 

and in the case B ± 

t x = et 2 = e£ 3 e 2r ? = he' 211 = sintf, t 5 = e = ±1, 
J = -et p = V 12 e 2r > = V 2l e- 2r > = costf, U = 2t p + sin ■& tan ti(e* - ee^) 2 , 

V u = V 22 = vf } = vf } = V 4 = 0, V 3 {1) = -VP = V 12 - V 2l . (5.47) 

with the free parameters rj. The parameter e = ±1 distinguishes between the models 
A , B , respectively. 

Let us introduce a replacement of the parameter rj that will become very useful in 
the following discussion. In the A case we shall set 

3 2„ el-^cosiT 1 



( (- COS 17 



while for the case B we use the definition 

4 — 1 COS 17 

U q (su(2) x su(2)) Symmetry. First of all we note that the AB Hamiltonian does not 
immediately match the expression of the general U g (su(2) xsu(2))-invariant Hamiltonian 
in Tab. 3. However, we can apply the transformation (5.12) which preserves integrability 
in order to restore the U g (su(2) x su(2)) invariance. 
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We find that the AB Hamiltonian can be brought into the form in Tab. 3 if some of 
the parameters are related: Explicitly we find that the parameters must be related 
by 

V n = V, V 3 {1) = -e 2ih qJ + e 2ih q-% - 2V, 

V 12 = e 2ih q~ x J + V, V 3 (2) = -e 2ih q- x J + e 2ih q~\ - 2V, 
V 21 = e 2lh qJ + V, V 3 {3) = -e 2i ^q~ x J + e 2lh qt p - 2V, 
V 22 = V V 3 W = -e 2if3 qJ + e 2lh qt v - 2V, 

V A = e 2ih (q + q- l )J - e 2lfl (q + q~ l )t p + AV, (5.50) 

with a new parameter V. Furthermore the tj, parameters must take the form 

(5.51) 



ti = to, t 2 = -e 2 ^H , t 3 = -e 2 ^q-H , U = e 2 ^qt , t 5 = -e 2 ^+ 2 ^, 



with a new parameter t . The parameters for the transformation take the form 

|o 1 = a 2 = -ie ai ^( g + 9 - 1 )- \U, 
\b 1 = b 2 = -\e 2i ^{q-q- 1 )-\U, 

a 3 = b 3 = 0, (5.52) 
and the twist parameters are constrained by 

Note that we have set ao = 1. We have also adjusted a 2 , b 2 such that H' — K' — 0. The 
resulting values for the parameters of Tab. 3 read 

A' = \e 2i ^(q + g-% + \U, D' = -\e 2i ^{q + g"% -\U + V, 

B > = _|e 2i /i(g + g-% + \U, E' = e 2i ^(q + q' 1 ).] + §e 2 ^(g + q~ l )t p -\U + V, 

C = -e ifl+if3 Kt Q (q + a" 1 ), F' = -e~ lh ~ ih K~H G (q + a" 1 ), 

Qi — e -ifi-Hh-if3 JJ — e +ifi+ih-ifa 

H' = 0, K' = 0. 

(5.54) 

The ten parameters A', . . . , V are thus given as functions of six parameters J, U, to, t p , V 
and k. Three further parameters are implicitly given by ao, a 2 (overall scaling and iden- 
tical shift) and b 2 (affecting H and K only) which preserve the form in Tab. 3 and which 
we have fixed above. One of the ten parameters cannot be chosen continuously due to 
the relation G'/L' = e 2ifl = ±1. 

Comparison of Hamiltonians. We find that in the cases A ± the symmetry can be 
restored only if e 2v = e or cos-$ = 0. Moreover, q — ±1. This means that we cannot 
explain the two-parametric model A ± in general using our methods. We shall therefore 
disregard this case here. In order to relate the cases A to our model we would have 
to find a more general transformation than (5.12) to recover the symmetry. It would be 
interesting to find out if this is possible. 
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Conversely, the models B 1 * 1 can be brought to the form in Tab. 3, and in what follows 
we shall exclusively consider these cases. We match the tk in (5.515.47) by demanding 

e = _ e *A+»* g = e V>/3. (5.55) 

The correct expressions for Vk in (5.50) follow by imposing the remaining four relations 
in (5.47) 

V = 0, J = -et p , £p + to = l, U = 2t p + (t p -t p 1 )(e(e 2r > + e - 2ri )-2). (5.56) 

Incidentally, these are precisely the constraints for our integrable Hamiltonian (5.8) when 
using the general AB coefficients (5.54). Therefore it should be possible to find suitable 
parameters x ± ,g in Tab. 4 to match the case B ± . 

We start the comparison by making the convenient choice 

A = -|7r(l + e), / 2 = / 3 = 0, e = -e 2if \ q = e 2 \ (5.57) 

To match G' = e - ^ 1 we have rescale our Hamiltonian in Tab. |3[4| 

H' hk = e -^H hk . (5.58) 

Furthermore we have to match our normalization C = F' by setting 

K = e~ ifl . (5.59) 
The ratio L'/G' = L/G = e 21 ? 1 = — e is determined through the x 



— i — s(x^) (30) 
G x~ — s(x~) 

This constraint is satisfied when (our) U 2 = — e or when q 2C = — e. Both values are 
possible, let us discuss them separately: 

The condition q 2C = —e holds for the following two pairs (x + ,x~) 

. iq 1 ± Jl -Aeg 2 (q 1 / 2 - eq- 1 / 2 ) 2 _ ei It a/ 1 - ^g 2 (q 1 ' 2 - eq- 1 / 2 ) 2 

fJ-> >^ Qft 

2g 1 - eq 2g 1 - eq 

(5.61) 

We match $ of B 1 * 1 using the coefficients A' — B' or C 

B' - A' e~ ifl B-A gie^-ee- 71 ) 2 

cos v = —et r - 



p e 2 V + e -2 V e 2r, + e ~2 V q ^1 - Aeg 2 - €e -v)2 ' 
smit = t = — — = — — — = =F • (5.62) 

e 2 V + e -2 V e 2r, + e ~2r ? Q T ^ _ 4^2^ _ V 7 

These two relations are compatible with the identity cos 2 d + sin 2 -d = 1 and they can be 
solved for g 

9 = 1 ■ (5.63) 

4 cosh(?7 + if i)y cosh 2 (77 + ifi) — cos 2 ■& 
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In the parametrization (5.49) the above relations simplify significantly 

q = e(i-ecos^)(e-cos^) l-gcosa x ± = iei ±i (5 64) 

9 cos^(e-e+^)(e-e-^)(e 2 -l)' q * £-costf ' x l< * " 

Conversely, the condition U 2 = — e holds for the following two pairs (x + ,x~) (note 
the change of sign w.r.t. (5.61)) 

+ _ iq 1 ± v 7 ! +4e^ 2 (g 1 /2 + eg -i/2)2 _ _ ez x ± ^1 + 4^2^1/2 + ^-1/2)2 
" 1 + eg ' ~ 2# 1 + eq 

(5.65) 

We match 1? of B using the coefficients A' — 5' or C" 

B'-A' e"^ 1 5 -A v/l-^ 2 (e 2 "-e-^) 2 

COS 77 = — 6t = = = zb — 

p e 2 ^ + e- 2? ? e 2 " + e- 2 " G y/l + 4ec/ 2 (e f ? + ee~") 2 ' 

sin # = t = — ~s s~ = — n 7^ = ± , • (5.66) 

e 2v + e- 2 '? e 2r > + e- 2? ? G ^/l + Aeg 2 (e^ + ee^) 2 

These two relations are compatible with the identity cos 2 1? + sin 2 — 1 and they can be 
solved for g 

9 = y • (5-67) 

4 sinh(r] + ifi)\ sinh 2 (?7 + ifi) + sin 2 1? 
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Again the parametrization (5.49) simplifies the above relations 

sini?(l — £ cosi9)(£ — cost?) 
costf(f - e+^KC - e-*>)((f + ^i) cost? - 2) ' 

1 — £ cos it 1 , ie£ sin it 1 _ ie sin 1? 

9 = e ^ "7 q" ' X = q c ' X = 1 1 q • 5 ' 68 

£ — COS 17 COS 17 — £ 1 — 4 COS 17 

Note that several different points in the constrained parameter space of x ,g corre- 
spond to the same Hamiltonian given in terms of d. Furthermore, we have not been very 
careful about selection of branches of U — VU 2 , q c = a/ q 2C in (2.62,2.64) as well as for 
the above square roots. We expect that these signs ambiguities are equivalent to shifting 
1? and irj by multiples of tt/2. 

Comparison of Bethe Equations. Here we will show how the 3^ models proposed by 
Alcaraz and Bariev |34| can be embedded into our U g (h)-spin chain model by comparing 
Bethe equations. 

Consider the Alcaraz-Bariev Bethe equations 34 

^ sinh(A fe - + ry) 

l _ » sinh(Aj 1 ) - Xj ± rj) f l sinh(A^ - - 2 V ) 
f = \ srnHX^ - A, - V ) f = \ sinh(Ai 1} - Xf + 2rj) ' 
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By inserting $(z) into the expression for exp(2Afc) in |34| we can write the relation 
between z k and A& as follows 

exp 2A fc ) = e~ 2ri -, — = HV^ t~ ■ 5.70 

The energy of an eigenstate is given by 

N 

E = YS*i + zJ 1 )- (5-71) 

3=1 

It is straightforward to match the auxiliary (second) Bethe equation with our Bethe 
equations (5.16) where fa = fa = by equating q = e 2r) and 

u(y k ) = exp(2A fc - 2A ) r , 

q-q 1 

w k = exp(2Ai 1} - 2A ) - y~ 1 ^^ ■ (5.72) 

q — q~ l 

Matching of the first Bethe equation requires 

z - c -ih a -c-i/2 lT -i yk-x + x+-e^q c+l / 2 Ux~z k 

We can bring the above relations (5.72.5.73) between X k and y k in the same form as 
(5.70) 

u(y k ) H = exp(2A ) vyks(y k ) (5.74) 

q-q 1 J q-q 1 

with 

_ 2 + z 2 k -e- i ^U- 1 q c +y 2 (2 + (U 2 -q- 2C )/q- 1 x+s{x+))z k + e- 2i ^q 
y k s(y k )-q x s(x ) ^ _ e _ i/lf/ „ lgC _ 1/2(f/2 - ? _ 2C) ^ - e - 2ihq -i ■ 

(5.75) 

We set q 2C = e 2 ^ 1 = — e and use (5.61). Noting that q 1 x + s(x + ) = qx s(x ) = — e the 
prefactors imply that we have to choose the shift parameter Ao as follows 

exp(2A ) = eig(q - q' 1 ) . (5.76) 

Some further manipulations then show that our relation (5. 74 5. 75) is precisely the same 
as the one from the Alcaraz-Bariev equations (5.70). Finally, the eigenstate energies E 
(5.71) match with our result E (5.19 5.17 5.23) using the transformation parameters 

a = e l/1 ^, ai = 2e^-, a 2 = e %h - , a 3 = 0. (5.77) 
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The Spectrum of the A 1 * 1 Cases. The Bethe equations derived by Alcaraz and 
Bariev are almost the same for the A and B^ 1 cases. Actually, the only difference is the 
definition of a global parameter r as a function of rj, namely 

cosh 2r = e sin 2 $ + cos 2 $ cosh 2i] vs. cosh 2r = cosh 1r\ (5.78) 

for cases A and B, respectively. This implies that set of solutions to the Bethe equations 
for model A with 7] = rj^ is the same as the set of solutions for model B with rj = rjB 
where tja and 7/b are related by 

e sin 2 $ + cos 2 cosh 2r]\ = cosh 2r] B (5.79) 

Moreover the dispersion relation is the same for both cases, and therefore the spectra 
of the models agree. Note that our reparametrization (5.48,5.49) of 77 has been chosen 
carefully so that the above relation holds for equal £ = £a = £b- 

In fact this leads to a puzzle concerning our above results: We have not been able to 
find a transformation of the Hamiltonian for model A that restores manifest U g (su(2) x 
5u(2)) symmetry. Therefore we have not been able to relate this Hamiltonian to ours. 
Nevertheless the matching of spectra implies that the two Hamiltonians are indeed related 
by a similarity transformation for a suitable choice of parameters. Consequently, there 
must exist a more elaborate transformation to make the U g (su(2) x su(2)) symmetry of 
the Hamiltonian of model A manifest and match it with our Hamiltonian. 



Comparison of Scattering Matrices. To shed more light onto the equivalence of 
the A^ and B 1 * 1 cases we shall compare the corresponding scattering matrices. These 
were derived in |34| for excitations c\ and c\ with momentum p = — ilogz above the 
(ferromagnetic) vacuum and they can be compared to the results of Sec. 4.2. We will use 
the parametrization (5.48,5.49) which allows us to compare the expressions directly with- 
out the further need to transform parameters. It turns out that the diagonal scattering 
matrix elements precisely match 

(S A )^i, z 2 ) = (S A )£(*i, z 2 ) = (S B )\l( Zl , z 2 ) = (S^tKz,, z 2 ). (5.80) 

The off-diagonal elements also agree up to a simple factor 



(S A )lj(z 1} z 2 ) = (S^lKz,, z 2 ) ^ , (SaKKz,, z 2 ) = (S B )H( Zl , z 2 ) . 



(5.81) 



with 



The fact that the off-diagonal elements differ by reciprocal factors of f{zi)/f(z 2 ) 
alone implies that the spectra are the same. Therefore there must exist a similarity 
transformation between the two Hamiltonians. In fact the transformation is simple for 
scattering states: Each excitation cj> with momentum p is multiplied by a factor of f(e ip ) 
while the other excitations c\ are left alone. It is not hard to convince oneself that this 
leads to the above relative factors between the S-matrix elements. However, it is not as 
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straightforward to express the similarity transformation for spin chain states because the 
the transformation must be non-local due to the momentum dependence of the factors 
f(e ip ). Astonishingly the Hamiltonian remains local after the similarity transformation 
which implies that the transformation must be of a very special kind. 



An Alternative Coalgebra? The above discussion has made it clear that the A ± 
and B 1 * 1 Hamiltonians are equivalent. Consequently the A Hamiltonian must obey a 
U 9 (su(2) x su(2)) symmetry albeit in a non-standard form. It is conceivable that this 
symmetry is generated by an alternative coproduct U q (i)). This is in fact a promising idea 
because it would explain the locality of the A^ Hamiltonian: Assuming that there exist 
an alternative (fundamental) R-matrix, we could right away derive a nearest neighbor 
Hamiltonian which is likely to be the one of the A^ cases. Thus it would be exciting 
to find an alternative coalgebra structure for U q (t)) and also find a suitable similarity 
transformation to the canonical one derived in this paper. Unfortunately this issue is 
outside the scope of the present paper. 



5.6 Quantum-Deformed EKS Model 

At the (x + ,x~) values (0,0), (00,00) as well as (s(0),s(0)), (s(oo), s(oo)) one obtains 
the quantum-deformation of the EKS model |25|. For instance at (x + ,x~) = (0,0) the 
parameters of the rescaled Hamiltonian Ti.'j k = (q 1 ^ 2 /G)Hj, k read 

A' = -B' = -D' = E' =l(q + q- 1 ), G = L~ l = q l/ \ C = F' = H' = K' = 0. 

(5.83) 

It is not too hard to see that after a twist this leads to the quantum-deformed EKS 
model, i.e. an standard integrable nearest-neighbor spin chain with U g (u(2|2)) symmetry 
and spins in the four-dimensional fundamental representation. 

Note that this model is somewhat singular. For instance the Bethe equations (5.16) 
do not apply directly, but the limit (x + ,x~) — > (0,0) has to be taken carefully. In 
particular, there are two allowed limits of the y k , either y k — ► or y k — > s(0). Let A x of 
them be of the first kind and A 3 = N — Ni of the second. They should scale like 

y k ~ exp(2X k 1) )qx~, y k+Nl ~ s(exp(2\ k 3) )qx~) . (5.84) 

Furthermore, the w k must scale like 

i q + q' 1 exp(2X k 2) )qx- 
2# q- q 1 9 2 {q - q 1 ) 2 

The limit of the Bethe equations (5.16) agrees with the standard form of U q (u(2|2)) 
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Bethe equations, cf. |46|, up to a twist 

l = e+m / sinh(Af-j) \ " ^ smh(A["-Af + ft) 
\smh(\M + ft) ) f.\ s inh(A<" - Af - ft) ' 
* Binh(Af - A'" + ft) " sinh(Af - Af - 2ft) ft s iih(A< 2 > - Af + ft) 
1 " l\ 8 inh(Af - Af - ft) Jj smh(A< 2 » - Af + 2ft) *i 8 mh(A<?> - A« - ft) ' 

" 3in h(Af-AfTft) 
1-6 M-MAf-Af-ft) <5 ' 86) 

with g = e 2? \ Here the two different kinds of have led to the additional Bethe equation 
needed for standard U g (u(2|2)) chains. 

Finally, the limit of the total momentum (5.17) and the total energy (5.19) yields 



e iP = e^s-JVi) TT smh(A^. 1) + h) 
f = \ sinh(Ai 1) - h) ' 



E = cosh(2h)K + Y sinh 2 (2ft) m . (5.87) 

jr[ sinh(A^ + h) sinh(A^ - h) 

We observe that only the Ajj. carry momentum and energy, the Ajj. do not contribute. 

5.7 Another Quantum-Deformation of the Hubbard Model 

Some curious values for the parameters (x + ,x~) appear to be (0,s(0)), (s(oo),oo) and 
(s(0),0) (oo,s(oo)). In the limit q — > 1 they obviously approach the values (0, oo) and 
(oo, 0) corresponding to the Hubbard model. One has to be careful in using these values 
because several common combinations of x ± turn out to be singular and have to be 
regularized properly, e.g. by taking a limit. 
For definiteness we shall use the pair 

(x + ,x-) = (0,s(0)); (5.88) 

the other combinations yield similar results. For convenience we rescale the Hamiltonian 

ia 1 / 2 

H' hk = ^- H j)k . (5.89) 

The regularized coefficients of this Hamiltonian read 

1 - 2(g - g- 1 ) V 



A' — B' — —D' = —E' 



Wi-(?-rW' 

a = F' = q + q-\ G' = +iq +l/2 , V = -iq' 1 ' 2 , H' = K' = 0. (5.90) 
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In particular, we find A' = B' and D' = E'. This is interesting because the pair-hopping 
terms of the oscillator Hamiltonian in Tab. 5 vanish just as for the Hubbard model (we 
set k = —i for convenience) 
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j,k 














+ iq 


H/2 t r (I _ 
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(5.91) 

Note that this Hamiltonian is not obviously hermitian. It would be interesting to find 
out if its spectrum is nevertheless real (for some non-hermitian twist). 

The closed chain with this Hamiltonian can be diagonalized by the Bethe equations 
(5.16). The regularized momentum relation (5.17) reads 

e.^V V-3-;+'^ , (5.92) 
1 - {q-q- Y ) 2 g 2 



and the energy of an eigenstate (5.19) takes the form 

E' = 1 ~ 2{<1 ~ q ~ l)2 jL= (K - 2N) + Y(iq +l / 2 e^ - iq-^e^). (5.93) 

42V 1 - (?- q~ l ) 2 g 2 ~t 

Note that this dispersion relation is not isotropic. This fact is related to our choice of 
x ± . E.g. choosing x such that the dispersion relation is isotropic leads to the Alcaraz- 
Bariev model discussed in Sec. 5.5, It may be desirable to investigate this model further 
to see if it has interesting (physical) properties. 



6 Conclusions and Outlook 

In this paper we have considered quantum deformations of the threefold central extension 
f) = psu(2|2) kM 3 of the Lie superalgebra psu(2|2). We have set up the Hopf algebra U g (h) 
together with its fundamental R-matrix and applied it to derive quantum deformations 
of the one- dimensional Hubbard model. 

A first important result is that the quantum deformation U g (l)) of the Lie superalgebra 
f) is possible despite its non-standard structure (central extension). Moreover, we have 
constructed an invariant fundamental R-matrix which obeys the Yang-Baxter equation, 
crossing symmetry and fusion (with a suitable choice of overall phase factor). This 
suggests that the Hopf algebra can be made quasi- co commutative and quasi-triangular. 
Curiously, the matrix structure of the fundamental R-matrix is determined by U,j(h) 
invariance alone. This feature is related to the representation theory Uq(h) which we 
have outlined and which is largely analogous to the undeformed superalgebras f) and 
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0u(2(2) . We have also discussed hermiticity conditions and we were forced to choose q to 
be real. Usually in quantum groups U 9 (g) one can also choose q from the unit circle (such 
that q + q" 1 is real), and it would be important to understand if and how hermiticity 
can be achieved in that case for our algebra. 

We have applied the nested Bethe ansatz to the R-matrix in order to diagonalize it 
and to write down the Bethe equations for a closed spin chain. The Bethe equations are 
structurally similar to the Lieb-Wu equations for the one-dimensional Hubbard model: 
One should view the former as the quantum deformation of the latter; the Lieb-Wu 
equations are of the algebraic (XXX-like) type, while ours are trigonometric (XXZ- 
like). A question for future work is whether an elliptic (XYZ-like) deformation of our 
fundamental R-matrix and the Lieb-Wu equations exists. 

In the final section we have derived an integrable Hamiltonian from the R-matrix 
and compared directly it to the one-dimensional Hubbard model and some of its gener- 
alizations. Our Hamiltonian has manifest U g (su(2) x su(2)) symmetry, three indepen- 
dent parameters and it can be further deformed in several canonical ways (rescaling, 
shifts, twists) while preserving integrability. By adjusting the parameters we were able 
to match the one-parametric Hubbard chain as well the cases B ± of a two-parametric 
Hamiltonian proposed by Alcaraz and Bariev |34|. We have also identified a potentially 
interesting combination of parameters for which our Hamiltonian has no pair-hopping 
terms. We have not managed to explain the cases A ± of the Alcaraz-Bariev Hamilto- 
nian in general, but we have provided an argument based on the Bethe equations why 
this should be possible nevertheless. It would therefore be important to show explicitly 
how to relate the cases A , perhaps by finding an alternative coalgebra for U g (f)(2)). 
Furthermore an investigation of the possible condensed matter theory applications of 
quantum-deformations of the one-dimensional Hubbard model is worth performing. Do 
these models also display signs of superconductivity like the undeformed model? It is 
also worth finding out if any of the other electronic models discussed in the introduction 
can be obtained from our Hamiltonian. 

Here we have investigated the Hopf algebra U 9 (fj), but the analysis is far from com- 
plete. In fact, the symmetry algebra in the undeformed setup is much larger, it appears 
to be a Yangian (double) |21|. The corresponding quantum deformation would be a 
quantum affine algebra. One will need this algebra in order to formulate the universal 
R-matrix of which the R-matrix derived in this paper is the fundamental representation. 
A classical and undeformed analysis |24| has shown that the affine algebra is not just 
fj[tt, m -1 ]. It is rather a deformation of u(2|2)[w, w -1 ] which includes only one tower of 
central charges and one tower of inner automorphisms. Therefore it is quite clear that 
this quantum affine algebra is not Kac-Moody, and one cannot directly apply the general 
framework associated with such algebras. 

Finally, we should ask how our results can be applied to the AdS/CFT correspon- 
dence. Quantum-deformed (XXZ-like) spin chains with U 9 (su(2)) and U 9 (su(3)) sym- 
metry have indeed appeared in this context, in particular for so-called beta-deformed 
M = 4 gauge theory |47| with complex deformation parameter 1 48-50 1. However, the 
argument used in |49| tells us that the symmetry considered in this paper cannot apply 
to conformal gauge theories such as |47|: The two bosonic subalgebras su(2) of psu(2|2) 
play different roles; one originates from the internal symmetry which is deformed, and 
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the other one originates from conformal symmetry which remains undeformed. In con- 
trast, quantum deformations typically apply to the whole of a Hopf algebra, as is the 
case for our algebra. Nevertheless, it seems possible to deform the internal su(4) sym- 
metry and the associated S 5 , so why should it not be possible to deform su(2, 2) and 
the associated AdS 5 7 Let us speculate about such a complete quantum deformation 
of the AdS/CFT correspondence: an "AdS g /CFT 9 " duality. The field theory would 
have to have quantum-deformed conformal and internal symmetries, and it would most 
probably be formulated on some sort of non-commutative spacetime. The string theory 
would be formulated on a quantum-deformed (AdS§ x S 5 ) q background, which would 
be both curved and non- commutative. Perhaps the most convenient definition would be 
in terms of coset spaces like SOg(6)/SO g (5) with non-commutative coordinates. Despite 
the expected dire technical complications, the planar limit of these dual models would 
stand a good chance of being integrable, and the world sheet S-matrix would then be 
given through our quantum-deformed R-matrix. It would be exciting to find out if this 
proposed picture is more than daydreaming. 
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A Some Useful Relations 

The fundamental representation is defined in terms of the variables x constrained by 
the quadratic relation (2.61) 

x + q _ 1 . . / x + qx~\ i 

— + \-qx +ig(q-q 1 )[ V =-• A.l 

q x + qx~ \qx~ x + J g 

An expression which depends on x can therefore be written in many ways which are 
equivalent upon (A.l). For example one may choose to eliminate all (x~) k with k ^ 
0, 1 (because the constraint is quadratic in x~) in order to get a unique representative 
within the equivalence class of the expression. Unfortunately this particular expression 
is typically not the most economical one. However such a representative can be used for 
the purpose to show that a certain expression is identically zero (or that two expressions 
are equivalent). 

In this appendix we would like to present various identities involving the x ± which 
were used in the paper or may prove useful otherwise. 
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One Set of x^. 



u = q 1 u(x + ) — — = qu(x ) + 
^9 



i q Tl x ± s(x ± ) 



2g g 



q-q- 



(A.2) 



q 1 x + s(x + )—qx s(x ) = —ig(q — q V )u+\{q + q 1 ) 

= -igq^(q - q^)u(x + ) + g" 1 
= ~igq(q ~ q~ l )u(x~) + g 

q c U - q- c \J- 1 
~ q-cjj _ q cij-i 



One y. 



y + s(y) y + s(y) 



One Set of x^ and a 



y- 



9 q-q~ 



x+ - s(y) s(x )-y _ 2( , , 



x - s(y) s(x+) - y 
x + — y s(x + ) — y u(x + ) — u(y) 



x~ — y s(x ) — y u(x~) — u(y) 



Two Sets of x^. 



X O SiX-\ 1 SiXry 1 X -\ S(Xr\ 1 SiX-\ 1 



-x 2 s(xj")-s(xj) x 1 - s(xt) s(xf)-x 2 



X 



X 2 X l s(x 2 ) _ od rr 2 tt r l ) ^(^2 ) 



^2~ - S 04 ) ' U(x 1 )-U(xt) 



(A.3) 



) = ( X+ ? + 1 g +S(X ) = 2C 2 

- s(x+) (x-) 2 + 1 + s(x+) g 1 ' J 

g/:r+ -g" 1 /^" = _i + 7 + x = _ 7 _x (A.5) 

3 C L r (^) 2 + 1 = q -Cfr-l (* + ) 2 + 1 = X+ ~ X ~ (A 6) 

q x- + s(x~) q x+ + s(x+) qCU-q-CU- 1 { } 



2/2 + 1 8(y)a + 1 =l-( 9 - 9 -W (A.7) 



«(y) = ±*^ (A.8) 



g 2G f/ 2 (A.9) 

(A.10) 



= q 2C ^U 2 1 U 2 2 (A.ll) 



i 2 ^ _ i = q 2Cl+2 U 2 x (A. 12) 



q^Ul^P^ ^( (A. 13) 
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(C 2l P^K' 2 ) (C^P^Ki) 




(C^P^K,) (C 2 ,P 2 ,K 2 ) 
Figure 4: Scattering process and transformation of central charges. 

B Braiding from Scattering Problem 

An equivalent approach to derive and incorporate the braiding factors for the coproduct 
(2.45) is based on using an S-matrix which acts as a permutation operator and commutes 
with the coproduct [A (5), S] = for each generator 3 G U 9 (f)). This relation is equivalent 
to cocommutativity (2.49) upon identifying the S-matrix with the R-matrix as usual 
S = VK. 

What we want to stress here is the fact that for the quantum deformed case the intro- 
duction of the braiding in the coproduct is necessary to allow for a quasi cocommutative 
Hopf algebra. If one would only take the quantum-deformed algebra with the coproduct 
(2.37) there cannot be any R-matrix which transforms A to A op for the central elements. 
We will identify the additional braiding element with the central charges in such way 
that the coproduct will indeed become cocommutative on the center. 

In particular, let us consider two short modules with central charges (Ci, Pi,Kx) and 
(C2, P21 K2). For the coproduct (2.37) we have 

A(«P) = <P <g> 1 + q 2C <g> qj, (B.l) 

but 

A OP OP) = 1 ® + <p ® g 2£ . (B.2) 

We need to change the coproduct by introducing an additional braiding factor to make 
it quasi-cocommutative. The importance of having a quasi-cocommutative Hopf algebra 
lies in the fact that they have an R-matrix which not only intertwines the modules, but 
satisfies the quasi-triangularity condition from which the Yang-Baxter equation follows. 

Pairwise Scattering. Now we will consider the representation structure when the S- 
matrix acts on chains of multiplets. First we consider the scattering matrix of two short 
multiplets (see Fig. 4|) 

S l2 : (Ci) ® (C 2 ) — > (C' 2 ) ® (C[). (B.3) 
For each generator 3 £ U 9 (f)) we want the S-matrix to be invariant 

[A(3),S] = 0. (B.4) 

In particular, this relation must hold for the central charges. 



50 



An obvious way to conserve the total charge C\ + C 2 is to demand that the individual 
charges are merely interchanged 

C[ = C U C' 2 = C 2 . (B.5) 

Taking the coproduct (2.37) of the other central charges ^3, K and evaluating them before 
and after scattering we have 

P 1 + P 2 q 2Cl = P2 + P[q 2C \ 

K iq - 2C ' 2 + K 2 = K 2 q~ 2Cl + K[. (B.6) 

Combining the short multiplet constraint (2.29) with (B.5) we have the following equa- 
tions to be solved 

P 1 K 1 = P[K[, P 2 K 2 = P' 2 K' 2 . (B.7) 

The substitution 

P 1 = 1 - q^Ut, K x = q- 2C ' - C/f 1 , 

P 2 = (1 - q^U 2 )U h K 2 = {q- 2C * - U 2 l )U^ (B.8) 
simplifies the calculations dramatically and we getp^j 

P[ = (1 - q^U x )U 2 , K[ = {q~ 2 ^ - U^)U 2 \ 

P^ = 1 - q 2C *U 2 , K' 2 = q- 2C * - U 2 \ (B.9) 

Factorized Scattering. A basic requirement for a factorized i^-particle S-matrix S n 
for any permutation 7r G Sk is that it forms a representation of the permutation group 
Sk, i-e. SnS^r = S-ktt 1 - It imposes certain relations on Pj and Ki. We consider the 
permutation that interchanges three modules as follows 

S w : (Ct) <g> (C 2 ) <g> ((%) >-> (4) ® ® (4)- (B.10) 

This process can be represented in two different ways. The factorized scatterings (23 — > 
3'2' and r'3' — > 3"1") lead to the following relations 

P^P 2 -rP^-P' 2 q 2G \ P[ = P, (B.11) 
for the process 23 — > 3'2' and 

= (Pi + P 2 g 2Cl ) - {P[ + P' 2 q 2C ^q 2C - A + P 3 g 2Cl+2C2 , 
P^' = P^ (B.12) 

for the process 1'3' -> 3"1". 



ln Thcre are two roots for the second order equation - one of them corresponds to trivial scattering 
and we did not mention it here. 
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Fusion. Furthermore, the overall process can be represented as a pairwise scattering 
of composite multiplet (C12) with (C3). The former multiplet is long but for particular 
values of C\ and C2 it splits into two short multiplets. The central charges P[ 2 3 and 
K[ 2 3 become related by the following constraints 

Pi + P 2 q 2Cl + P 3 q 2Cl+2C ' 2 = P3 + P[q 2Ca + P^ 2Cl+2Cs , 
Kig -2C 2 -2C 3 + R2q -2C S +Rs = K ^ q -2 Cl -2C 2 + ^/ g -2C 2 + ^/ (R13) 

The results (B.12) and solution of (B.13) should correspond and it constrains the form 
of central charges. An educated guess consists in choosing them in the following way 

fc-i 

P k = ga(l-q 2C *U 2 k )l[U 2 , 

i=i 

k—l 

K k = %-^-U^)Hur 2 . (B.14) 

a 

which completely agrees with (2.52). 

The following remark is noteworthy here. We have derived the dependence of the 
central charges P and K on the braiding factors U by demanding the condition that 
scattering with a pair of particles can be treated as a successive scattering with one 
particle and then with another one. Alternatively we can treat U as a function of P and 
K and put UU^ 1 = 1 as a constraint (which is for sure satisfied in (B.8) automatically). 
If this constraint is not satisfied then the Hopf algebra for f) = su(2|2) x IR 2 cannot be 
quasi-triangular. We can drop one central charge (say K and thus have su(2|2) x M. 1 
algebra) then this relation does not constrain us any longer and the Hopf algebra based 
on su(2|2) K M 1 is quasi-triangular. 

Fundamental S-Matrix. We write down the explicit form of the fundamental S- 
matrix S = VIZ in Tab. 6, Its coefficients are the same as for the fundamental R-matrix 
given in Tab. 2. 
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